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ABSTRACT 


An  approximate  quasi-three-dimensional  method  of  solution  is  de¬ 
scribed  for  the  sudden  area  expansion  flow  process  in  the  relative 
coordinate  system  under  the  assumptions  of  frictionless,  adiabatic  flow 
of  a  perfect  gas.  The  objective  was  to  determine  the  radial  variation 
of  the  flow  properties  at  the  trailing  edge  of  a  rotor  given  the  measured 
flow  properties  downstream  of  the  rotor  through  use  of  the  streamtube 
approximation.  Results  are  derived  from  one  configuration  of  the  blunt 
trailing-edge  supersonic  compressor  rotors  tested  at  AEDC.  The  re¬ 
sults  were  felt  to  be  an  unsatisfactory  representation  of  the  average 
flow  conditions  at  the  trailing  edge,  and  it  is  shown  that  the  reasons  for 
this  are  relatable  to  neglecting  the  free  turbulent  shear  flows  occurring 
in  the  blade  wakes  in  conjunction  with  neglecting  the  large  radial  second¬ 
ary  flows  which  apparently  occur  in  the  flow  field  of  the  rotor  investi¬ 
gated.  It  is  shown  that  the  streamtube'  approximation  can  produce 
grossly  inaccurate  results  when  free  turbulent  shear  is  neglected  on  the 
streamtube  boundaries.  A  proposal  for  continued  work  is  also  given. 
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SECTION  I 
INTRODUCTION 


The  description  of  the  three-dimensional  flow  in  turbomachines  is 
so  mathematically  complex  and  unrewarding  to  date  that  analysis  of  test 
data  relating  to  turbomachines  is  often  accomplished  by  assuming  a 
model  of  the  flow  through  a  blade  row  and  solving  the  design  equations 
at  stations  between  blade  rows  where  there  are  no  blade  forces  and 
where  the  analysis  may  be  carried  out  in  the  absolute,  or  stationary, 
coordinate  system  (see  Chapter  III  of  Ref.  1,  for  example).  As  a  model 
approximating  the  influences  of  a  blade  row,  the  analyst  may  choose 
pertinent  cascade  results  and  apply  empirical  correction  factors  to 
account  for  the  stacking  position  of  the  blade  element  in  the  blade  span 
and  the  additional  frictional  effects  occurring  in  the  blade  end  regions. 
Alternately,  he  may  apply  a  completely  mathematical  model  that  satis¬ 
fies  the  equations  of  fluid  flow  in  two  dimensions  or  one  that  basically 
originated  from  results  of  cascade  testing  or  testing  of  rotors  of  simi¬ 
lar  design  and  then  apply  the  additional  three-dimensional  correction 
factors. 

Once  the  analyst  has  chosen  his  model,  his  results  are  completely 
relatable  to  that  model.  If  the  results  agree  with  the  experiment,  then 
the  model  is  considered  good  for  that  application,  but  the  model  may  be 
limited  only  to  that  application.  If  the  results  only  partially  agree  with 
experiment,  the  disagreement  is  attributable  to  influences  which  the 
analyst  chose  to  neglect  or  could  not  account  for. 

This  type  of  approach  has  its  value  to  compressor  or  turbine  de¬ 
signers  although  it  may  leave  the  real  flow  process  in  the  rotor  unex¬ 
plained.  At  AEDC,  under  the  sponsorship  of  ARL,  several  configura¬ 
tions  of  a  high-reaction  shock-in-rotor  compressor  rotor  have  been 
tested  (Refs.  2  and  3).  The  blade  element  is  of  blunt  trailing-edge  design 
for  the  purpose  of  promoting  supersonic  constant-area  diffusion  to  sub¬ 
sonic  trailing-edge  velocities  while  eliminating  the  starting  difficulties 
generally  associated  with  earlier  designs  of  supersonic  compressor 
blading  (Ref.  4). 

The  process  of  supersonic  constant-area  diffusion  cannot  yet  be 
handled  by  purely  theoretical  means.  However,  the  process  can  be 
approximated  because  experimental  results,  such  as  those  in  Refs.  5 
and  6,  show  that,  if  the  passage  length-to-height  ratio  is  optimum,  then 
the  losses  in  the  passage  approach  the  losses  of  a  single  normal  shock 
at  an  equivalent  passage  inlet  Mach  number.  Thus,  at  least  for  optimum 
design,  the  trailing-edge  flow  conditions  can  be  approximated. 
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Between  the  rotor  trailing  edge  and  the  measuring  station  down¬ 
stream  of  the  rotor,  unusually  large  mixing  losses  were  anticipated 
because  of  the  large  trailing-edge  thickness  of  the  blading.  The  quasi- 
two -dimensional  studies  in  Ref.  7  showed,  however,  that  if  the  passage 
flow  process  behaved  properly,  then  the  Mach  number  at  the  trailing 
edge  would  be  small  enough  so  that  the  additional  loss  due  to  the  abrupt 
increase  in  flow  area  would  be  of  acceptable  magnitude. 

The  results  for  the  rotors  tested  at  AEDC  have  been  far  below  that 
expected  of  the  blading,  and  the  author  has  attempted  in  Ref.  8  to  de¬ 
termine  the  cause  of  the  poor  performance  through  re-evaluation  of  the 
basic  flow  model.  It  was  determined  that  the  passage  diffusion  process 
was  incomplete  due  either  to  the  blade  passage  design  or  to  the  occur¬ 
rence  of  separation  or  both  of  these.  Therefore,  the  trailing-edge  flow 
property  distribution  was  uncertain  since  the  basic  passage  flow  model 
could  not  be  effectively  applied,  and  the  real  influence  of  the  abrupt 
area  increase  at  the  trailing  edge  was  unknown. 

In  an  attempt  to  learn  more  about  the  flow  property  distribution  at 
the  trailing  edge,  and  thereby  more  about  the  flow  process  in  the  rotor, 
it  was  decided  to  apply  the  flow  model  in  a  slightly  unusual  way.  Rather 
than  proceeding  with  the  flow  from  upstream  to  downstream  of  the  rotor 
as  attempted  in  Refs.  9  and  10,  the  measured  conditions  upstream  of  the 
rotor  were  used  to  estimate  shock  loss  alone,  whereas  the  measured 
conditions  downstream  of  the  rotor  were  used  to  calculate  the  conditions 
at  the  trailing-edge  plane.  This  report  relates  the  consequences  of 
using  the  measured  downstream  conditions  in  such  a  manner. 

Initially,  it  was  suspected  that  the  quasi-two-dimensional  theory 
(Ref.  7)  would  be  limited,  if  not  unsatisfactory,  because  considerable 
radial  mass  transfer  was  evident  in  the  results  of  testing.  Therefore, 
the  quasi-three -dimensional  technique  presented  in  this  report  was 
developed  both  to  allow  radial  shifting  of  stream  surfaces  in  the  flow 
field,  a  feature  not  allowed  in  the  much  simpler  quasi-two-dimensional 
theory,  and  to  allow  determination  of  the  radial  variation  of  flow  proper¬ 
ties  at  the  trailing  edge. 

Although  the  method  presented  in  this  report  does  permit  these 
goals  to  be  attained,  the  results  must  be  considered  quite  unsatisfactory 
as  a  representation  of  the  average  flow  properties  at  the  trailing-edge 
plane.  The  method  involves  solution  of  the  integral  equations  of  motion 
in  a  relative,  or  rotating,  coordinate  system.  Flow  is  assumed  friction¬ 
less,  and  it  is  shown  that  this  assumption  is  the  primary  reason  for  the 
unsatisfactory  results.  However,  in  contrast  to  general  occurrences, 
shear  stress  is  shown  to  be  the  determining  factor  throughout  the  flow 
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field,  not  just  in  the  region  of  hub  and  tip  casing  walls.  It  is  shown  that 
the  frictionless  equations  of  motion  will  predict  extremely  large  losses 
in  relative  total  pressure  when  used  in  the  manner  described  in  this  re¬ 
port,  and  it  is  further  shown  that  free -turbulence,  coupled  with  large 
radial  secondary  flows,  is  a  reasonable  explanation  for  the  unsatisfac¬ 
tory  results. 


SECTION  II 

FORMULATION  OF  THE  EQUATIONS  FOR  SUDDEN  AREA  EXPANSION 
WITH  RADIAL  VARIATION  OF  FLOW  PROPERTIES 


The  bLunt  trailing-edge  blading  may  be  expected  to  have  substantial 
additional  losses  in  relative  total  pressure  due  to  the  viscous  mixing 
process  which  occurs  after  the  abrupt  area  increase  at  the  trailing  edge. 
The  analysis  in  Ref.  7  develops  quasi-two-dimensional  equations  to 
estimate  the  additional  loss.  The  results  indicate  that  the  losses  due 
to  the  sudden  area  expansion  process  are  relatively  small  if  the  axial 
Mach  number  at  the  trailing-edge  plane  is  not  too  dose  to  unity.  It  has 
been  indicated  in  Section  I  that  the  passage  of  the  blunt  trailing-edge 
compressor  was  designed  as  a  constant-area  diffuser  to  reduce  super¬ 
sonic  inlet  relative  velocities  to  subsonic  velocities  at  the  trailing  edge. 
Therefore,  minimum  sudden  expansion  loss  depends  on  the  proper 
design  of  the  passage  and  requires  a  compromise  between  longer  pass¬ 
age  length  to  produce  lower  exit  Mach  numbers,  thus  increasing  friction 
loss  in  the  passage,  and  shorter  passages  to  reduce  machine  weight, 
thereby  increasing  sudden  expansion  loss. 

The  quasi-two-dimensional  analysis  for  the  loss  in  relative  total 
pressure  due  to  the  sudden  area  expansion  process  cannot  determine 
the  radial  variation  of  loss;  it  can  only  give  average  values  for  the 
annulus.  Furthermore,  the  quasi-two-dimensional  theory  requires  that 
the  radial  velocity  be  zero  for  the  entire  process,  but  the  test  results 
have  indicated  that  there  is  considerable  radial  change  of  stream  surfaces 
through  the  wheel,  especially  at  design  speed.  The  use  of  the  quasi-two- 
dimensional  theory  implies  that  all  radial  shift  must  occur  within  the 
rotor. 

These  considerations  place  strong  limitations  on  the  use  of  the 
quasi-two-dimensional  theory;  and  a  method  has  been  derived  which 
should  determine  the  radial  variation  of  important  flow  properties  and 
the  approximate  position  of  streamlines  at  the  axial  location  of  the  blade 
trailing  edges.  The  method  uses  the  integral  flow  equations  for  inviscid 
flow  and  provides  solutions  for  the  flow  properties  as  a  linear  function 
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of  radius  over  an  arbitrary  increment  of  radial  annulus  height.  The 
method  allows  a  tangential  component  of  velocity  to  exist,  but  the 
solution  does  not  determine  tangential  flow  property  variations;  i.  e. , 
flow  is  still  assumed  axisymmetric  so  that  the  solutions  are  quasi- 
jhree- dimensional  in  nature,  representing  tangential  averages  at  any 
radial  station. 

The  method  employs  the  integral  equations  of  motion  referred  to  a 
relative  coordinate  system  rotating  with  the  angular  velocity  of  the  com¬ 
pressor  rotor.  As  discussed  in  Ref.  11,  the  choice  of  the  relative 
coordinate  system  leads  to  equations  involving  volume  integrals  of  the 
force  due  to  Coriolis  and  centripetal  accelerations.  As  stated  in  Ref.  11, 
these  volume  integrals  generally  cannot  be  evaluated  without  complete 
knowledge  of  the  flow  properties  in  the  entire  control  volume.  In  some 
cases,  flow  property  measurements  have  been  made  that  allow  an 
approximate  calculation  of  the  volume  integrals.  These  cases  are  used 
to  determine  weighting  factors  applied  to  the  surface  integrals  at  en¬ 
trance  to  and  exit  from  the  control  volume  such  that  the  volume  integrals 
determined  by  use  of  the  weighting  factors  are  equivalent  to  the  volume 
integrals  of  the  approximate  calculations. 

The  use  of  the  equations  of  motion  in  a  rotating  coordinate  system 
should  generally  be  limited  to  a  region  where  every  fluid  element  is 
participating  in  curvilinear  motion  imparted  by  the  angular  velocity  of 
the  wheel.  In  turbomachinery,  its  use  should  be  limited  to  flow  inside 
a  rotating  member  where  the  force  exerted  by  the  blades  generally  domi¬ 
nates  the  frictional  effects  at  the  nonrotating  surfaces  that  tend  to  change 
Jhe  transportational  velocity  of  a  fluid  element.  In  the  present  application, 
j;he  equations  of  motion  for  frictionless  flow  in  a  relative  coordinate  sys¬ 
tem  are  being  applied  to  a  region  downstream  of  the  rotor,  but  friction  on 
the  nonrotating  surfaces  of  the  casing  walls  could  exert  considerable  in¬ 
fluence  on  the  solution  because  the  flow  property  measurements  are  the 
result  of  a  frictional  flow  field.  In  the  present  solution,  frictional  effects 
are  ignored,  and  it  is  assumed  that  the  inertial  forces  maintain  the  trans¬ 
portational  motion  of  every  fluid  element. 

It  is,  of  course,  possible  to  view  the  flow  field  in  any  coordinate  sys¬ 
tem  that  may  be  chosen.  The  relative  coordinate  system  was  chosen  be¬ 
cause  flow  is  more  likely  to  be  steady  after  the  rotor  in  such  a  system.  ' 

In  the  nonrotating,  or  absolute  coordinate  system,  the  flow  may  not  be 
considered  steady  after  the  rotor.  Thus,  unsteady  effects,  as  well  as 
frictional  effects,  would  have  to  be  considered  for  complete  solution  of 
the  problem  in  the  absolute  coordinate  system. 
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2.1  DESCRIPTION  OF  THE  PROBLEM 

Consider  a  rotor  of  blunt  trailing-edge  blades  and  locate  a  rotating 
cylindrical  coordinate  system  with  its  origin  on  the  axis  of  the  rotor  and 
its  axial  unit  vector  (k)  pointed  along  the  axis.  The  rotational  velocity 
of  the  coordinate  system  (&)  will  be  assumed  to  match  the  rotational 
velocity  of  the  rotor.  A  positive  value  of  the  axial  unit  vector  (k)  points 
in  the  positive  direction  of  Si  in  the  right-hand  sense.  A  positive  value 
of  the  tangential  unit  vector  (afg)  points  in  the  direction  of  rotation. 
Defining  the  radial  unit  vector  (sfr)  to  point  away  from  the  axis  when  it 
is  positive  sets  up  the  right-hand  cylindrical  coordinate  system  as  shown 
in  Fig.  1. 


At  a  particular  distance  from  the  axis  (r)  a  cylinder  of  radius  r  is 
made  to  intersect  the  blades.  If  this  cylinder  is  then  spread  onto  a 
plane,  the  blade  profiles  may  be  represented  as  shown  in  Fig.  2.  This 
is  the  general  representation  of  compressor  blades  in  cascade.  The 
blades  are  pictured  with  a  spacing  s,  chord  length  c,  and  trailing-edge 
thickness  th. 

A  plane  perpendicular  to  the  axis  containing  the  blade  leading  edges 
in  Fig.  1  is  represented  by  the  line  z  =  Z2  in  Fig.  2.  A  similar  plane 
containing  the  blade  trailing  edges  is  represented  by  z  =  Z3  in  Fig.  2. 

The  plane  z  =  Z4  in  Fig.  2  is  located  far  enough  downstream  so  that 
complete  mixing  of  the  flow  has  occurred. 

Figure  3  represents  a  plane  containing  the  axis  of  the  rotor  inter¬ 
secting  the  bounding  walls  of  the  rotor  casing.  This  is  the  general 
representation  of  the  compressor  in  the  meridional  plane.  The  curve 
C-C'  represents  the  generatrix  of  the  compressor  casing.  The  distance 
from  the  axis  to  the  compressor  casing  (r^)  may  be  an  arbirtary  function 
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pf  z,  the  axial  distance  from  the  origin  of  the  coordinate  system.  The 
purve  D-D'  represents  the  generatrix  of  the  compressor  hub,  and  its 
fiistance  from  the  axis  (r^)  may  be  an  arbitrary  function  of  z. 


Fig.  2  The  Cascade  Representation  of  the  Blading 


Fig.  3  A  View  of  the  Rotor  in  the  Meriodional  Plane 
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The  region  (R)  of  interest  for  the  solution  of  the  losses  due  to  sudden 
area  expansion  will  now  be  defined.  The  region  (R)  is  enclosed  by  the 
planes  z  =  Z3  and  z  =  Z4  and  the  surfaces  generated  by  C-C'  and  D-D' 
between  these  planes.  The  region  will  have  a  total  surface  area  {A). 

Fluid  flow  may  enter  the  region  only  through  the  surfaces  area  A3, 
where  Ag  is  defined  as  the  annulus  area  at  z  =  z3  excluding  the  area  A3' 
occupied  by  the  blade  trailing  edges.  Fluid  flow  may  leave  the  region 
only  through  the  annulus  area  A4  at  z  =  Z4. 

It  is  assumed  that  the  flow  through  area  A3  at  z  =  Z3  has  not  yet  ex¬ 
perienced  the  sudden  area  expansion  and,  also  that  it  is  in  such  a  posi¬ 
tion  that  the  blades  may  exert  only  an  infinitesimal  moment  on  the  fluid. 

The  problem  is  to  define  the  equations  necessary  to  determine  the 
relationship  between  the  fluid  properties  at  the  planes  z  =  Z3  and  z  =  Z4. 

Since  the  problem  involves  a  great  many  unknowns  and  since  the 
radial  variation  of  some  of  the  flow  variables  cannot  be  specified  with 
certainty  at  the  plane  z  =  Z3,  it  is  advantageous  to  divide  the  region  (R) 
into  a  selected  number  of  streamtubes.  An  example  of  one  such  stream- 
tube  is  shown  in  Fig.  3  bounded  by  the  generatrices  A- A'  and  B-B\ 

Then,  with  some  confidence,  radial  variation  of  certain  flow  variables 
within  each  streamtube  may  be  ignored.  As  an  illustration,  the  adiabatic 
energy  equation  requires  that  the  total  enthalpy  be  constant  along  a  stream¬ 
line.  In  some  cases,  the  assumption  of  constant  total  enthalpy  within  the 
entire  region  may  be  crude.  However,  if  the  region  is  divided  so  that 
radial  variation  of  total  enthalpy  is  small  within  each  streamtube,  accuracy 
may  be  improved. 

The  streamtube  approximation  is  inaccurate  whenever  there  is  net 
mass  transfer  between  streamtubes;  i.  e. ,  secondary  flows  within  the 
region.  This  represents  a  limitation  in  the  solution  of  the  problem. 


2.2  ASSUMPTIONS  FOR  THE  SIMPLIFICATION  OF  THE  EQUATIONS 

In  order  to  simplify  the  equations,  the  following  assumptions  were 
made: 

1.  The  fluid  is  a  perfect  gas.  This  assumption  allows  the 
use  of  the  equation  of  state  for  a  perfect  gas, 

p  =  pRt  (1) 


7 


AEDC-TR-71-1 1 1 


where  p  is  the  static  pressure,  p  is  the  density,  t  is 
the  static  temperature,  and  R  is  the  gas  constant. 
Furthermore,  the  specific  heat  at  constant  pressure 
<Cp)  is  constant  so  that  the  ratio  of  specific  heats  {«) 
is  a  constant  which,  thereby,  tends  to  simplify  com¬ 
putation. 

2.  The  weight  of  the  fluid  in  the  region  (R)  is  negligible. 

In  comparison  to  the  inertial  forces,  the  weight  of  the 
gas  is  usually  quite  small,  and  this  assumption  allows 
a  volume  integral  to  be  neglected. 

3.  Flow  is  steady.  It  is  necessary  to  consider  a  rotating 
coordinate  system  for  rotor  flows  in  order  to  efficiently 
use  this  approximation.  The  integral  equations  to  be 
applied  to  this  problem  contain  volume  integrals  of  un¬ 
steady  effects  requiring  complete  knowledge  of  the 
flow  at  every  point  in  the  region  at  every  instant  of 
time.  It  is  usually  not  experimentally  possible  to 
evaluate  this  integral  for  the  type  of  problem  under 
consideration.  Time  averaged  values  of  the  flow  prop¬ 
erties  reduce  the  error  in  this  assumption. 

4.  Flow  is  adiabatic  and  frictionless.  The  energy  equa¬ 
tion  will  normally  include  the  heat  generated  by  friction 
along  the  path  of  a  fluid  particle  if  the  flow  is  assumed 
steady  and  if  heat  transfer  to  the  boundaries  of  the 
region  is  ignored  when  it  is  assumed  that  the  heat  due 
to  friction  is  immediately  and  entirely  dispersed 
throughout  the  fluid  element.  However,  the  mathemat¬ 
ical  complexity  of  including  the  frictional  term  in  the 
momentum  equatipn  requires  the  total  neglection  of 
friction  for  the  present  calculations. 

5.  FLow  is  axisymmetric  throughout  the  region  (R).  This 
assumption  implies  that  there  is  no  variation  of  the  flow 
properties  in  the  tangential  (stg)  direction  at  any  partic¬ 
ular  axial  position  (z)  and  radial  distance  from  the 
axis  (r).  If  the  plane  z  =  is  far  enough  downstream, 
this  assumption  produces  very  little  inaccuracy  at  this 
station;  however,  because  of  the  unusually  large  thick¬ 
ness  of  the  blade  trailing  edges,,  this  assumption  becomes 
progressively  poorer  as  the  plane  z  =  Z3  is  approached. 

The  plane  at  z  =  Z3  is  essentially  made  up  of  two  different 
types  of  surfaces  that  partially  bound  the  region  (R).  One 
is  the  effective  surface  of  the  blade  trailing  edges  com¬ 
prised  of  the  actual  trailing-edge  surface  area  and  the 
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portion  of  the  blade  surface  boundary  layer  that  is 
intersected  by  the  plane  z  =  Z3.  The  area  of  this 
surface  was  termed  earlier  A3'.  By  assumption,  no 
flow  may  enter  the  region  through  this  surface,  but 
secondary  flows,  perhaps  caused  by  centrifugal 
forces  for  instance,  could  conceivably  cause  compo¬ 
nents  of  the  velocity  to  exist.  The  other  type  of  sur¬ 
face  at  Z3  comprises  the  area  between  the  blades 
through  which  flow  enters  the  region  (R),  and  it  was 
earlier  termed  A3.  In  general,  sizable  tangential 
variation  of  all  flow  properties  may  be  expected  in 
the  z3  plane,  particularly  if  the  blades  are  cambered, 
although  one  might  expect  these  variations  to  be  peri¬ 
odic  from  blade  to  blade  if  the  flow  relative  to  the 
blades  is  steady.  The  equations  used  to  describe  the 
sudden  expansion  process  will  account  for  the  assump¬ 
tion  that  the  velocity  is  identically  zero  over  the  sur¬ 
face  A3  ,  whereas  the  assumption  of  axial  symmetry 
will  be  used  to  neglect  any  tangential  variation  over 
the  surface  A3.  These  assumptions  lead  to  equations 
that  still  contain  terms  which  allow  possible  differ¬ 
ences  between  the  static  pressures  on  the  surface  A3' 
and  A3. 

The  results  of  the  calculation  presented  in  this  report 
are  derived  from  equations  which  assume  that  the 
static  pressure  on  the  blade  trailing-edge  surface  and 
the  static  pressure  in  the  free  stream  are  identical. 
However,  some  reports  (for  example.  Ref.  12)  indi¬ 
cate  that  there  may  be  considerable  difference  between 
these  pressures,  and  the  proper  equations  accounting 
for  this  are  presented.  In  this  case,  a  relationship  of 
the  form 


P3'  =  P3  f  (M3z»  A4/A3,  k ) 

is  required  where  P3'  is  the  pressure  on  the  trailing- 
edge  surfaces,  P3  is  the  free-stream  static  pressure, 
and  f(M3z,  A4/ A3,  * )  may  be  an  empirically  determined 
function  of  the  free-stream  axial  Mach  number  {M3Z), 
the  flow  area  ratio  (A4/A3),  and  the  ratio  of  specific 
heats  (*). 

6.  Secondary  flows  are  neglected.  Such  an  assumption  is 
likely  to  be  unrealistic  in  fluid  flow  as  complicated  as 
that  in  turbomachines .  By  way  of  example,  consider 


(2) 


9 


AEDC-TR-71-111 


the  space  behind  a  blade  trailing  edge.  In  all  like¬ 
lihood,  such  a  region  contains  an  unusually  large 
amount  of  low  energy  boundary  layer  material  com¬ 
pared  with  blading  of  more  conventional  design. 

This  low  energy  fluid  can  be  mbre  easily  centri¬ 
fuged  outward  to  .collect  at  the  casing  than  the 
main  stream  flow.  This  particular  type  of  sec¬ 
ondary  flow  is  probably  inherent  in  the  overall 
picture  of  flow  with  the  blunt  trailing-edge  blading. 

Various  other  types  of  secondary  flows  are  de¬ 
scribed  in  Chapter  XV  of  Ref.  1.  Experimental 
evidence  related  in  Ref.  1  indicates  that  wheels  of 
high  hub-to-tip  radius  ratios  may  be  expected  to 
reduce  secondary  flows.  In  any  case,  only  quali¬ 
tative  analysis  from  experimental  data  for  sec¬ 
ondary  flows  is  available  to  date  because  of  the 
mathematical  complexity  of  the  problem. 


2.3  THE  EQUATIONS  FOR  A  ROTATING  COORDINATE  SYSTEM 

The  following  is  a  list  of  the  equations  applied  for  solution  of  the 
sudden  area  expansion  phenomenon.  For  the  most  part,  they  are  taken 
directly  from  Ref.  11.  In  review,  the  derivation  of  these  basic  equa¬ 
tions  hypothesized  a  flow  that  is  representable  by  field  properties  that 
are  continuous  and  continuously  differentiable.  The  angular  velocity  of 
rotation  has  been  assumed  constant;  separation  at  the  planes  of  entrance 
to  and  exit  from  the  control  volume  is  assumed  not  to  occur;  and  flow 
may  not  both  enter  and  leave  the  control  volume  through  the  same  open 
surface  area. 

2.3.1  Momentum  Equation 


The  integral  momentum  equation  for  steady  relative  flow,  neglect¬ 
ing  friction  and  the  weight  of  the  gas,  is 

/W4  dm,\4  -  J  W3dmj4g  =  J  -npdA  -  J  pfl  x  (2W  +  H  x  r)  dv  (3) 

(A4)  (A3)  (A)  (R) 

where  dm^.  is  the  absolute  value  of  the  increment  of  the  mass  flow  rate 
defined  by 

dm.\i  =  (4) 
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where  p  is  the  density,  n  the  unit  vector  normal  to  the  element  of  sur¬ 
face  area  dA^,  and  Wj  is  the  relative  velocity  given  in  cylindrical  coordi¬ 
nates  by 

=  arWr[  +  (5) 

In  Eq.  (3),  the  areas  A3  and  A4  are  the  areas  through  which  the  fluid, 
respectively,  enters  and  leaves  the  region  (R).  The  region  (R)  is 
bounded  by  a  surface  of  total  area  (A).  By  convention,  the  unit  normal 
vector  (n)  points  away  from  the  region  (R)  so  that  the  negative  sign  in 
the  integral  of  the  static  pressure  (p)  indicates  that  the  normal  stress 
is  exerted  on  the  surface  element  dA  by  the  surroundings.  In  general, 
the  final  volume  integral  of  Eq.  (3)  for  the  Coriolis  acceleration 
(2Sz  x  W)  and  the  centripetal  acceleration  (fl  x  (il  x  ?))  cannot  be  evalu¬ 
ated  accurately  because  of  lack  of  sufficient  flow  measurements  within 
the  region. 

2.3.2  Continuity  Equation 


The  continuity  equation  may  be  written  in  integral  form  for  steady 
flow  in  the  relative  coordinate  system  as 

J  dm  a  3  =  /dm  A4 

*  A  3 )  (A  4)  (6) 


where  dm^.  is  given  in  Eq.  (4). 


2.3.3  Energy  Equation 


The  energy  equation  for  steady,  adiabatic  flow  in*a  rotating  coordi¬ 
nate  system  is 


d'RH'  =  0 


where  d'p(  )  represents  the  change  along  the  path  of  a  fluid  particle  in 

the  rotating  coordinate  system  and  where  H'  is  the  relative  total  enthalpy. 
For  a  perfect  gas,  this  leads  to 


Ti  =  r; 


K  -  1 

2kR 


n2(r23  -  r») 


The  streamtube  approximation  is  used  in  approximating  an  average 
value  of  Tg'  in  each  streamtube  based  on  the  mass  averaged  value  of  T4' 

the  streamtube  center  at  Z4,  and  an  arithemetic  average  of  the  inner  and 
outer  boundary  radii  at  Z3. 


11 


AEDC-TR-71-1 1 1 


2.3.4  First  Moment  of  Momentum 


The  axial  component  of  the  first  moment  of  momentum  for  steady, 
frictionless  flow  of  a  weightless  gas  in  a  rotating  coordinate  system  is 
given  by 


jjl'z  =  /  k  .  r3  x  W3dmA3 
(A3) 


J  k  •  r4  x  W4  dmA<  -  f  k  ■  r  x  [ptt  x  (2W  +  fi  x  r)]  dv 
(a4) 

(7) 


If  no  force  is  exerted  in  the  tangential  direction,  this  equation  equals 
zero. 


2.4  DERIVATION 


The  problem  which  now  remains  is  the  rearrangement  of  the 
above  equations  into  the  most  useful  form.  To  do  this,  several  addi- 
tional^relationships  are  necessary.  The  relative  Mach  number  vec¬ 
tor  (M')  is  given  by 


where 


-*  re  . 

M  =  — 

a 


a  =  \'KRt 


Psing  the  isentropic  relation  (Ref.  11)  gives 


-  =  1  +  -Hi  M2 

t  2 


The  static  temperature  may  be  replaced  by  the  relative  total  temperature 
(T')  when  the  Mach  number  in  the  above  relation  is  considered  the  rela¬ 
tive  Mach  number  (M').  Then  the  local  speed  of  sound  is  given  by 

yJ~K\\ T7 


f\  +  M'2'* 


I'f 


so  that  the  relative  velocity  (W)  may  be  replaced  by 

W  - 


M  V'kR  T' 


(,  +  !^L  .'*)’■ 


(8) 


It  is  now  convenient  to  define  the  geometry  of  the  surfaces  bound¬ 
ing  the  region  (R).  The  effective  thickness  of  a  blade  at  the  trailing 
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edge  (th)  may,  in  general,  be  a  function  of  the  radial  distance  from  the 
axis.  If  th  is  considered  the  arc  length  at  a  particular  radius  r,  then 
the  angle  intercepted  by  this  thickness  is  th/r.  The  cross-sectional  area 
at  the  trailing  edge  occupied  by  a  single  blade  is  then 

r3t  th/r 

A  =  J  J  rd0dr 

f  3h  0 

r3 1 

=  J  th  dr 
'3h 

where  the  subscript  3  indicates  the  plane  at  z  =  z3,  the  plane  of  the 
trailing  edge,  t  implies  the  tip  or  outer  casing,  and  h  implies  the  hub 
or  inner  casing.  Then  for  n  blades  the  incremental  area  at  the  trail¬ 
ing  edge  (dA3')  is  given  by 

dA3 '  =  n  rh  dr  (9) 

For  the  assumption  of  axisymmetric  flow,  the  element  of  area  at  z  =  z3 
through  which  fluid  may  pass  (dAg)  is  given  by 

dA3  =  (2nr  -  n  lh)dr  (10) 

The  element  of  surface  area  A4  at  z  =  Z4  is 

dA*  =  2rrrdr  (ID 


for  axisymmetric  flow. 

To  the  present  point  in  the  derivation,  the  control  volume  has  been 
considered  that  volume  bounded  by  the  axisymmetric  surfaces  of  the 
inner  and  outer  casing  walls,  the  plane  of  the  rotor  blade  trailing  edges 
(z3),  and  the  downstream  measuring  plane  (Z4).  The  streamtube 
approach  is,  next  applied  to  divide  this  control  volume  into  smaller 
regions.  This  procedure  is  necessary  to  allow  the  linear  representa¬ 
tions  of  the  flow  properties  in  each  streamtube  to  be  sufficiently  accurate 
and  to  reduce  inaccuracy  in  the  integration  of  the  energy  equation.  At 
the  downstream  measuring  plane,  several  stream  surface  locations  may 
be  determined  by  use  of  the  measurements  at  this  station  in  conjunction 
with  the  continuity  equation.  It  is  assumed  in  the  present  application 
that  the  hypothetical  stream  surfaces  of  each  streamtube  are  parts  of 
right  circular  cones  between  measuring  planes  and  the  plane  of  the  blade 
trailing  edges.  A  significant  portion  of  the  problem  involves  locating 
the  stream  surfaces  at  the  blade  trailing  edge  plane.  It  is  assumed  that 
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there  is  no  net  mass  transfer  across  the  stream  surfaces.  Two  exam¬ 
ples  of  such  stream  surface  configurations  are  shown  in  Fig.  4  as  a 
cross-sectional  view  containing  the  rotor  axis.  Figure  4  (top)  shows 
the  boundaries  of  a  streamtube  that  might  occur  when  there  are  no 
extra  measurements  at  an  axial  plane  between  the  rotor  trailing  edge 
plane  and  the  downstream  measuring  plane.  Figure  4  (bottom)  shows  a 
possible  configuration  when  the  extra  measuring  plane  is  included  for 
determination  of  the  weight  numbers  for  the  volume  integrals  (see 
Section  2.8). 
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Fig.  4  Examples  of  Arbitrary  Stream  Surface  Configurations 

In  order  to  define  the  elements  of  surface  area  for  the  axisym- 
metric  stream  surface  bounding  the  streamtubes,  it  should  be  recalled 
that  curves  like  A- A',  B-B',  C-C'  and  D-D'  in  Fig.  3  represent  the 
generatrices  of  axisymmetric  surfaces.  A  streamtube  like  that 
bounded  by  A-A'  and  B-B"  in  Fig.  3  has  one  surface  boundary  (B-B  ) 
that  is  closer  to  the  axis  of  generation  which  will  be  termed  the  inner 
boundary  (subscript  I)  and  one  surface  boundary  (A-A')  that  is  farther 
away  from  the  axis  of  generation  which  will  be  termed  the  outer  bound¬ 
ary  (subscript  O).  The  elements  of  area  for  the  inner  and  outer  stream- 
tube  boundaries  are  given  by,  respectively, 

dAi  =  27rr[(z)  dsj 
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and 


dAo  =  2?rro(z)  dso 


where,  in  general,  rj  and  tq  may  be  arbitrary  functions  of  z  and  ds  is 
the  element  of  length  between  z  =  Z3  and  z  =  Z4  in  a  plane  containing  the 
axis.  For  a  straight  line,  the  length  elements  ds  may  be  given  by 


ds| 


dso 


<1  z 
COS 

d  y 

ens  £0 


where  ?  is  the  angle  of  a  tangent  to  the  surface  measured  with  respect 
to  the  axis  in  a  plane  containing  the  axis.  The  elements  of  area  of  the 
the  streamtube  surface  boundaries  then  become 


dAt 


2)7i [(/)  dz 
cos 


and 


dAo  - 


2~  *()(./-)  Hz 
fus  fo 


(12) 


(13) 


The  unit  normal  vectors  to  each  element  of  the  total  surface 
area  (A)  may  now  be  determined.  First,  the  planes  z  =  Z3  and  z  =  Z4 
are  perpendicular  to  the  axis  of  the  rotor;  therefore,  the  normals  to 
these  surfaces  are  parallel  to  the  unit  vector  k.  The  unit  normal  vec¬ 
tors  point  away  from  the  region  (R)  so  that  the  unit  normal  vector  at 
z  =  Z3  is  equivalent  to  -k.  Similarly  at  z  =  Z4,  the  unit  normal  vector 
equals  +k.  To  designate  the  components  of  the  vectors  normal  to  dAj 
and  dAo  consider  Fig.  3.  The  sign  of  the  angle  §  which  the  wall  makes 
with  positive  axial  direction  is  positive  if  r4  >  r3  for  Z4  >  Z3  and  nega¬ 
tive  if  r4  <  r3  for  Z4  >  Z3.  Figure  5a  shows  two  possible  local  orienta¬ 
tions  of  the  inner  boundary  generatrix.  The  normal  to  the  inner  stream- 
tube  boundary  is  given  in  both  cases  by 

n  -  k  sin  £\  -  ar  cos  (14) 

Similarly,  Fig.  5b  shows  two  possible  local  orientations  of  the  outer 
boundary  generatrix.  The  normal  to  the  outer  streamtube  boundary  is 
given  in  both  cases  by 


n 


ar  cos  -  k  sin 


(15) 
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The  element  of  volume  for  the  region  (R)  is  given  by 

dV  =  -27TT  dr  dz  (16) 

for  axisymmetric  flow  where  r  may  be  a  function  of  z. 


The  various  vector  products  occurring  in  Eqs.  (  3)  and  (  7)  may  be 

evaluated  using  Eq.  (5),  the  relation 

-»  -*  -» 
r  =  arr  +  kz 

and  the  fact  that  ft  =  kft.  Then  the  following  relations  result: 


O  x  W  =  -arQW(9  -  aQ  ftWr  i 

ft  x  (ft  x  r)  =  -ar  ft2r  ■ 

£  •  r  x  ^  ( 

k  -  r  x  [/>ft  x  (2W  +  ft  x  r)]  =  2p(2rtyr  / 


(17) 


a.  Inner  Boundary 


b.  Outer  Boundary 


Fig.  5 


Geometric  Relation  of  the  Vectors  Normal  to  the  Axisymmetric  Stream  Surfaces 


The  following  relations  are  derived  from  Fig.  6  and  allow  substi¬ 
tution  for  the  components  of  the  relative  velocity. 


W'z  =  W  cos  y  cos  J3 

Wr  =  IV  sin  y 

W#  =  -W  cos  y  sin  $ 


(18) 
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Fig.  6  The  Three-Dimensional  Velocity  Triangle 

By  using  Eqs.  (1),  (4),  and  (5)  and  relations  in  Eqs.  (8)  to  (18),  it  is 
possible  to  rearrange  Eqs.  (3),  (6),  and  (7)  to  the  following  form 

Radial  Component  of  the  Momentum  Equation 


r30 

J  P3  Ms'2  sin  y3 

'31 


cos  y3  cos 


n  th  \ 
2m  ) 


rdr  =  A'  +  B'  +  Rf  +  R2(19) 


where 


r40 

A'=  J*  P4  M42  sin  y4  cos  y4  cos  /S4  rdr 
'41 


(20) 


1  14 

B'  =  —  /  (porO  -  Pin)  dz 

*3 


(21) 


R 


1 


n 

rrv'xR 


— — —  M'2  )  cos  y  sin  dv 
2  / 


Ra 


O2 

2jtkR 


(22) 


(23) 
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Tangential  Component  of  the  Momentum  Equation 


r3° 

J  P3  M32  cos2  y3  sin  (83  cos  ^83  ^1  -  rdr  =  C'  +  R3 

'31 


where 


r40 


/2 

P4M '4  cos2  yi  sin/34  cos  /S4  rdr 


r4! 


R-J  ^ 


-  0 


(R)  ^  V  2  J 


Axial  Component  of  the  Momentum  Equation 


r3°  .  .  . 

/  P3  (Hj  s  cos2  yt  cos2  ft,  +  f  )  (1  -  If)  rdr  =  D'  +  E'  +  F' 


r3 1 


where 


r40 

D'  h  J*  p4  (M  42  cos2y4  cos2  £4  +  -~Vdr 


r4I 


E'  H 


F'  = 


—  J  (piri  tan  £1  -  poro  tan  £0)  dz 
*3 

r30 

1  f  n  th  , 

J  Ps  rdr 


r3T 


Continuity  Equation 


(24) 


(25) 


(26) 


(27) 


(28) 

(29) 

(30) 


Vff  /  PaM,'co.y,co.fl,(l  +  ^Ms'1)  (l  -  f))  ■*  '  -|j  (31) 

r3I 


where 


r4° 

G  '  =  J  P4M4'  cos  y4  cos  (84  (l  -  —  M  32^  ^  rdr 


(32) 


f4t 
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Axial  Component  of  the  Law  of  Moment  of  Momentum 


J  p3  M  32  cos2  y3  sin  /S3  cos  £3(1  -  r2dr  =  II '  +  84  (33) 

r3 1 


where 


r*o 

H  -  /  p*m;2  cos^  y4  sin  /34  cos  /34  r2  dr  (34) 

r4 1 

R4'  h  _^Lr  /  _£(l  +  Hi  M'2)^  sinyrdv  (35) 

Equations  (19),  (24),  (27),  (31),  and  (33)  describe  the  process  of 
sudden  area  expansion  for  a  rotor  of  blunt  trailing- edge  blading  under 
the  present  assumptions.  In  order  to  actually  solve  the  equations,  it  is 
necessary  to  make  further  assumptions  and  to  select  five  particular 
quantities  to  be  obtained  by  the  solution.  The  following  sections  deal 
with  these  problems. 


2.5  FURTHER  ASSUMPTIONS  TO  SIMPLIFY  THE  PROBLEM 
2.5.1  Assumption  1 


The  quantities  A'  of  Eq.  (20),  C'  of  Eq.  (25),  D  of  Eq.  (28),  O’  of 
Eq.  (32),  and  H*  of  Eq.  (34)  involve  only  the  properties  of  the  flow  at  the 
plane  z  =  Z4.  These  quantities  may  be  evaluated  numerically  from  meas¬ 
urements  at  the  Z4  plane. 

For  the  specific  case  of  the  present  analysis,  it  is  assumed  that 
the  radial  component  of  Mach  number  is  identically  zero  at  z  =  z^; 
i.  e. ,  74  =  0  (see  Eq.  (18)).  Measurements  were  obtained  under  this 
assumption,  and  in  Ref.  8,  it  was  demonstrated  that  neglecting  radial 
velocity  produces  only  small  errors  at  this  plane. 

2.5.2  Assumption  2 

The  quantities  B'  of  Eq.  (21)  and  E'  of  Eq.  (29)  involve  the  pres¬ 
sure  on  the  streamtube  bounding  surfaces  and  the  shape  of  these  sur¬ 
faces  between  zg  and  Z4.  In  general,  neither  of  these  is  fully  known, 
and  therefore,  they  must  be  approximated. 
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Static  pressure  is  measured  only  on  the  casing  walls;  but  with  the 
assumptions  that  static  pressure  varies  linearly  in  the  radial  direction 
and  that  the  stream  surfaces  form  portions  of  right  circular  cones,  it 
is  possible  to  separate  the  functions  into  integrals  that  remain  constant 
for  all  streamtubes  and  coefficients  of  these  integrals  that  are  deter¬ 
mined  by  a  streamtube's  location  and  geometry.  The  following  discus¬ 
sion  shows  the  exact  procedure  leading  to  this  simplification. 


It  will  be  assumed  that  the  surfaces  of  a  streamtube  form  a  por¬ 
tion  of  a  right  cone  which  has  an  axis  coinciding  with  the  axis  of  the 
rotor.  This  means  that,  in  the  representation  of  Fig.  3,  the  stream- 
tube  has  straight  lines  denoting  stream  surface  shape  between  planes 
z  =  Z3  and  z  =  Z4.  Then  radial  distance  of  a  stream  surface  at  any 
axial  location  z,  Z3  <  z  <  Z4,  is  given  by 

7.\  —  Z 

rj(z)  =  r<ij  -  (r4i  -  r3j)  -  ,  i  =  I  or  0  (36) 

H  -  7.3 


The  slope  in  the  z  direction  of  a  stream  surface  is  approximated  by 

r4i  **  r3i 

tan  fi  =  - ,  i  =  I  or  0  (37) 

Z4  -  Z3  V  ' 


The  static  pressure  is  usually  measured  at  static  taps  positioned 
along  the  solid  boundaries.  The  variation  of  static  pressure  from  inner 
to  outer  walls  is  assumed  linear.  Then  for  a  particular  axial  station  (z), 
the  variation  of  static  pressure  with  radius  is  given  by 


p(r,  z)  =  P]w(z)  4  [pi(z) 


Ph(z)] 


r  ~rhU) 
rt(z)  -  rh(z) 


(38) 


where  r^(z)  and  r^-(z)  are,  respectively,  the  hub  and  tip  casing  wall 
radii  at  the  particular  axial  station  z,  and  where  p^fe)  and  pt(z)  are  the 
static  pressures  measured  at  the  hub  and  tip  walls,  respectively,  for 
the  particular  axial  station  z.  Since  it  is  of  particular  interest  to  know 
the  static  pressure  on  the  stream  surface,  the  streamtube  surface  radius 
of  Eq.  (36)  must  be  substituted  for  r  in  Eq.  (38).  Then  the  static  pres¬ 
sure  on  the  streamtube  surface  i  (=  O  or  I)  at  the  particular  axial  posi¬ 
tion  z  is  given  by 


p;(z)  =  Ph(z)  +  lpt(z) 


rh(z) 

rh(z)_ 


(39) 


-  tptU)  -  Ph(z)l 


z4  -  z  1  f  r4i  -  r3i 
*4  ~  *3  J  Lri(z)  -  rh(7). 
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Substitution  of  Eqs.  (36)  and  (39)  into  Eq.  (21)  leads  to  the  follow¬ 
ing  representation  of  B': 


B'  =  O40  “  r4l)  B{  +  (r^o  -  r4])  B2  +  2  (r3or40  ~  r3Ir4l)B3 
+  (r3o  -  rai)  B4  +  (r|o  -  r||)  B5 

where 


(40) 

(41) 

(42) 

(43) 

(44) 

(45) 


Similarly,  substitution  of  Eqs.  (36),  (37),  and  (39)  into  Eq.  (29)  yields 


E'  =  (rj,  -  r340) 


B'2 


z4  -  ^3 


+  030  ~  r|i ) 


B4' 


*4  “  z3 

+  (r30r 40  -  r3Ir'U ) 


+  (r 2.  -  ^  *  CrSo  -  £>  7^ 

z4  -  Z3  z4  z 

+  (r30r40  ~  r3Ir4j) 


(46) 


4  “  z3 


2  _ 2  \  F-2 _  ^  I _ _  _  _ ^  E3 


z4  “  z3 


+  (r30r4  0  “  r3lr4l) 


z  4  “  z3 


where  B \$  B'2,  B'4,  and  B'5  are  given,  respectively,  by  Eqs.  (41), 
(42),  (44),  and  (45)  and  where 


E9  = 


(47) 

(48) 


Pt  ~  Ph 
rt  -  Th 


(49) 
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It  should  be  noted  that  the  integrals  represented  by  B  i  =  1,  . . . , 

5,  and  E'^,  i  =  1,  2,  3,  involve  only  the  measured  values  of  the  wall 
static  pressure  and  the  annulus  wall  radii.  Therefore,  these  quantities 
may  be  integrated  without  knowledge  of  the  stream  tube  configurations 
and  considered  constants  in  obtaining  the  solution.  However,  since 
evaluation  of  the  integrals  requires  knowledge  of  the  annulus  wall  static 
pressure  at  the  plane  of  the  trailing  edge,  a  particular  problem  arises  at 
hub  wall  where  no  wall  static  pressure  may  be  measured  in  practice.  In 
the  present  application,  the  other  measured  static  pressures  along  the  hub 
wall  are  used  to  extrapolate  for  the  missing  value.  It  is  felt  that  since  the 
integrals  are  based  on  several  actual  measured  values,  little  error  is  in¬ 
volved  in  this  necessary  approximation. 


2.5.3  Assumption  3 


The  volume  integrals  R’l  of  Eq.  (22),  R'2  of  Eq.  (23),  R3  of 
Eq.  (26),  and  R'4  of  Eq.  (35)  require,  in  general,  complete  knowledge  of 
the  flow  at  every  point  within  the  region.  However,  it  will  be  assumed 
that  the  value  of  the  integrals  may  be  approximated  by  a  weighted  average 
of  the  values  of  the  integrands  at  the  planes  z  =  Z3  and  z  =  Z4.  This  is  a 
strenuous  approximation  since  it  completely  neglects  the  flow  process 
between  these  planes.  However,  the  weights  are  determined  from  solu¬ 
tions  of  the  equations  in  each  streamtube  when  the  surface  integrals  are 
known  at  the  measuring  plane  between  the  rotor  trailing  edge  and  the 
downstream  measuring  plane.  The  method  of  determining  the  weights 
will  be  discussed  in  Section  2.8. 


Under  this  assumption  and  the  assumption  of  axial  symmetry,  the 
quantities  R'i,  E'2,  R'3,  and  R  4  may  be  rewritten  in  the  following  form: 


20  (z4  -  z3)  \  1  -  W  -f30 


-LZjLi.  fr3°  P3M3COS  y,  sin  )3,(l  4-  —  Mj*Vl  -  ,  !?>£( 

Vrj  Jr3)  \  2  /\  2rrr  /  y/rj  ) 


(50) 


where  X  is  the  integral  evaluated  at  z  =  Z4 

r4l 


4  /  /  \  ’ 

X'  s  J  P4M4  cos  y\  sin  £4  fl  +  ~  ~  1 — M4  2  J  rdr 

r40 


(51) 


and  where  is  the  weight  associated  with  X 


R'2  - 


fi2  (z3  -  Z4)  j  ]  -  vv 2 


kR 


f-?  r  -( ■  •■■)  (■  •  •£■)  *■  •  ¥1 

'  r3 1 


(52) 
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where  Y  is  the  integral  evaluated  at  z  =  Z4 


r40  1 

Y'-J  P4  (] 


1  ^  m;3 


^  r^dr 


(53) 


>■41 


and  where  W2  is  the  weight  associated  with  Y'. 

20  (zj  -  z4)  (1  -  W,  r30 


Ri  - 


(1-W3  r30  /  k-  1  j  n  th  \  Vis'/. '  / 

-== -  \ — f  P3M3  sin  >'3  (1  +  -  M32)  (] - )  rdr  +  — 

VKH  }  VT3  JSI  V  2  /  V  2~r  )  v  1  i  S 


(54) 


where  Z'  is  the  integral  evaluated  at  z  =  Z4 


/*°  /  V  x 

=  J  P4M4  sin  M  +  — - —  M^2)  rdr 

r4|  \  2  / 


(55) 


and  where  Wg  is  the  weight  associated  with  Z'. 


20(z3-74)  (i-w*  /30  ...  /  K-l 

R4  =  - ;== -  ■)—_==-  f  p3M3  siny3  h - 

Vkr  I  vi  3  Jr3l  V  2 


M 


(56) 


where  Z"  is  the  integral  evaluated  at  z  =  Z4 

7"  =  J  P4M4  sin  y4  4  -  2  1  M  ;2^ 


r2dr 


(57) 


r4I 


and  where  is  the  weight  associated  with  Z". 

2.5.4  Assumption  4 


The  quantity  F  'of  Eq.  (30)  could  be  calculated  by  an  iteration  pro¬ 
cess.  According  to  Eq.  (2),  P3'is  dependent  on  the  axial  Mach  number, 
the  area  ratio,  the  free -stream  static  pressure,  and  the  ratio  of  specific 
heats.  Initially,  it  is  assumed  that  pg'  and  the  free-stream  static  pres¬ 
sure  (pg)  are  equal.  Solution  of  the  equations  under  this  assumption 
yields  the  information  for  the  evaluation  of  P3'by  Eq.  (2).  This  is  then 
used  to  calculate  F',  and  solution  of  the  equations  is  again  obtained. 

The  value  of  pg'is  then  recalculated  according  to  Eq.  (2)  and  the  pro¬ 
cess  continues  until  the  values  of  P3'for  two  successive  solutions  differ 
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by  less  than  some  convergence  criteria.  To  adapt  this  procedure, 
Eq.  (27)  is  better  written  in  the  following  form: 


/30  I”.," 

J  P3  M3' 

r3 1  L 


cos2  y3  cos2  J83 


,  r3  0 

rdr"K^riJ  P3thdr 

r3I 


*  t] 

1  n  .r3° 

=  D'  +  F/  -  K  —  —  f  p3'  til  dr 
k  2n  J 

r3I 


(27') 


where  for  the  initial  solution  K  =  0.  0  to  indicate  the  terms  involving  p3 
and  p3'  cancel  and  for  each  successive  solution  K  =  1.  0.  The  results 
presented  in  this  report  are  derived  for  K  =  0.0. 


By  using  Eq.  (40)  for  B',  Eq.  (46)  for  E',  Eq.  (50)  for  R'i, 

Eq.  (52)for  'R'2,  Eq.  (54)  for  R'3,  and  Eq.  (56)  for  R'4,  it  is  possible 
to  rewrite  Eqs.  (19),  (24),  (27),  and  (33)  into  a  form  for  which  the 
right  hand  side  is  independent  of  any  unknown  parameters  at  plane 

z  =  z3- 

f  P3M32  sin  }'3  cos  y3  cos  /33  ^1  -  — — ^  rdr  +  2(r3ir4|  -  r3or4o)  ^3 
r3i 

-  (r3I  -  r3o)  B  i  -  (r3I  -  r3o  )  B  5 
20(23-74)  .'3  0 

(1  -  Wi) 

r3 1 


VkRTj 


2.)  -'30  /  k-1  \)4  /  nth\ 

=  (1  -  W 1)  J  P3VI3  cos  >'3  sin  /S3^l  +  — - — *  M  32  J  ^  ^  JT& 


*  (1  .  lt2,  *  iUi  »;*)  (1  -  1±) 

kRT3  *3I  V  2  /  \  2nr  / 

-A'  +  (r40  ~  r4l)  Bj  +  (r40  -  rft )  B2  H  — *3-  WiX'  +  -  —  W2Y  '  (19') 

Vkrt;  krt  4 


'3  0 


r3I 


J  P3M3  2  cos2  y3  sin  $3  cos  £3 

20(74  -  73)  -30 


rdr 


1(84  ~  * 3)  (1  -  *3)  f  P3M3'  sin-y3(j  +  Ma2']  ’  (l  -  — ]  rdr 

VkRTs  '31  '  2  '  '  2trr  ' 


20(7  3  -  7.1) 

C'  l-  — — ■■■  —  W3Z  ' 
VkRT4 


(24  ) 
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/3o  r 

J  P3  ^3 

r3[  L 


2  cos2  y3  cos2  03 


Bl 


rdr  +  (r3Sf_rf0)  — +  (rfl  -  rf0)  - *— 

z4  “  z3  z4  “  z3 


+  (r3ir4|  -  1-3^40) 


z4  “  z  3 


2  1  +  (T3ir4I  “  r30r40 )  — — —  +  (r3ir4I  “  r3or4o)  3 


z  4  ~  z  3 


z 4  “  z3 


-  K  - — — 

k  2  n 


+  (r24i  -  r40) 


-r30 

J  P3  th  dr  =  D '  +  (r43  -  t^q) 


BZ 


r3I 
®  1 


z4  ~  z3 


r,  1  n  rri0 

—  K  —  —  J  P3 '  th  dr 


z4  -  z3  k  2 n 


(27") 


r3I 


r30  /  ,  v 

J  P3M3'2  cos2y3  sin  £3  cos  03  P  “  r2<Jr 

r3I  ' 

2fi(z4  -  z3>  (]  _  w<)  J  30  p3  sin  y3  ^  *  111  JJ|2y 


1  - 


V  kR  T  3 


n  th  \  „ 

- \  r2dr 

2rrr  / 


ir  +  w4z" 

VkRT4 


(33') 


Note  in  Eq.  (27”)  that,  although  p3'  is  not  actually  a  known  quantity,  the 
iteration  procedure  described  previously  would  ensure  a  value  of  pg'  for 
each  iteration  step. 

Equations  (19'),  (24  ),  (27”),  (31),  and  (33')  represent  the  final 
form  of  the  five  equations  describing  the  sudden  area  expansion  pro¬ 
cess.  The  following  section  will  describe  the  particular  assumptions 
necessary  to  limit  the  unknowns  on  the  left  hand  side  to  five. 

2.6  ASSUMPTIONS  CONCERNING  THE  PARTICULAR  FUNCTIONAL 
FORM  OF  THE  FLOW  VARIABLES 


2.6.1  Assumption  1 

It  has  been  assumed  that  static  pressure  varies  linearly  with  radius 
across  the  annulus  in  order  to  evaluate  integrals  of  streamtube  boundary 
pressure  forces.  The  flow  variables  are  to  be  considered  linear  func¬ 
tions  of  radius  in  each  streamtube  so  that  radial  changes  from  stream- 
tube  to  streamtube  may  be  represented  by  a  continuous,  piece  wise -linear 
curve.  These  two  approximations  are  consistent  only  if  the  static 
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pressure  {P3)  is  assumed  linear  across  the  annulus  and  not  piecewise- 
linear;  therefore,  at  z  =  Z3  the  static  pressure  at  any  radius  is  deter¬ 
mined  from 

P3  =  Cl  +  C2r,  ri,{z3)  <  r  <  ri  (z3)  (58) 

where 

c,  „ 

rt  (z3>  “  rK  (*3) 


and 


Cl  -  pt  (z3)  -  C2  rt  (z3) 


The  wall  static  pressure  is  measured  at  the  tip,  and,  in  order  to 
specify  the  variation,  it  must  be  determined  at  one  other  point,  such 
as  the  hub  wall.  Such  a  value  has  been  extrapolated  from  other  meas¬ 
urements  on  the  hub  wall  to  evaluate  the  integrals  of  the  resultant  force 
due  to  pressure  on  the  stream  surfaces,  and  this  pressure  must  be 
used  to  yield  the  necessary  function  in  the  solutions  here.  However, 
this  is  a  very  strenuous  assumption  since  the  extrapolated  pressure 
does  not  necessarily  obey  any  of  the  laws  of  fluid  flow.  The  extrapo¬ 
lation  is  felt  to  yield  reasonable  but,  certainly,  unjustifiable  estimates 
of  the  hub  wall  static  pressure  at  the  plane  of  the  rotor  trailing  edge. 

By  considering  the  simplification  that  this  approximation  provides  in 
allowing  direct  integration  of  B'j,  i  =  1,  2,  . .  .  ,  5,  and  B\,  i  -  1,  2,  3, 
considering  the  limited  number  of  equations  available  to  determine  the 
flow  variables  at  z  =  Z3,  and  considering  the  other  simplifications  neces¬ 
sary  to  this  problem,  then  the  extrapolation  to  obtain  p^(z3>  is  justified; 
but  it  is  noted  that  improvement  in  this  assumption  should  be  considered 
necessary  in  a  more  advanced  solution  of  this  problem. 

2.6.2  Assumption  2 

The  relative  flow  angle  at  the  blade  trailing  edge  03)  is  assumed  to 
vary  linearly  in  each  streamtube;  i.  e. , 

£3  =  b|  +  b2r,  r3j  <  r  <  r3o  (59) 

Use  of  the  integral  equations  that  are  being  applied  to  the  sudden  expan¬ 
sion  flow  process  demands  that  there  exist  no  flow  separation  at  either 
the  entrance  to  or  the  exit  from  the  bounded  flow  region  between  Z3 
and  Z4.  Therefore,  this  demand  is  obeyed  if  the  solution  for  j3 3  deviates 
very  little  from  the  average  angle  of  the  blade  trailing  edge  O3'). 
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Furthermore,  cascade  tests  have  shown  that  £3  is  greater  than  £3'  ex¬ 
cept  in  end  wall  regions  where  secondary  flows  may  cause  £3  to  be  less 
than  £3'.  These  considerations  are  expected  to  form  one  basis  for  judg¬ 
ing  the  applicability  of  the  present  method  for  solution  of  the  sudden  ex¬ 
pansion  flow  process. 

2.6.3  Assumption  3 

The  sine  of  the  relative  pitch  angle  (73)  is  assumed  to  vary  linearly 
in  each  streamtube;  i.  e. , 

sin  Y3  =  gl  +  g2r>  T3I  <  r  <  r30  (60) 

From  the  assumption  in  derivation  that  flow  may  only  enter  the  region 
at  z  =  Zg,  the  value  of  sin  7g  is  limited  to  values  between  +1  and  -1; 
i.  e.,  -90°  <  73  <+90°. 

2.6.4  Assumption  4 


The  relative  Mach  number  (M  g) 
each  streamtube;  i.  e. , 


is  assumed  to  vary  linearly  within 


M3  =  dj  +  dgr. 


i"3I  <  r  <  no 


(61) 


2.6.5  Assumption  5 

The  thickness -to-radius  ratio  (th/r)  may  be  considered  the  angular 
distance  per  blade  at  a  fixed  radius  through  which  no  flow  may  enter  the 
control  volume.  In  real  flows,  this  angular  distance  would  be  composed 
of  the  angle  subtended  by  both  the  blade  trailing  edge  and  some  flow  de¬ 
fect  that  is  due  to  boundary  layer  accumulation  on  the  blades.'  However, 
the  governing  equations  are  frictionless  and  axisymmetric  and  contain 
no  mechanism  that  would  determine  this  angular  distance.  The  equations 
certainly  do  prescribe  a  certain  flow  area,  but  they  give  no  preference  as 
to  its  linear  dimensions,  whether  tangential  or  radial.  This  implies  that 
either  th/r  or  the  stream  surface  locations  at  the  blade  trailing-edge 
plane  must  be  specified.  In  the  present  solution,  th/r  is  specified  as  the 
blade  trailing-edge  thickness -to-radius  ratio,  while  the  stream  surface 
locations  are  demanded  from  the  solution. 

Since  a  certain  flow  area  is  required  by  the  equations  and  since  the 
measurements  used  to  calculate  various  terms  in  the  equations  are  re¬ 
sults  of  real  flows,  it  is  reasonable  to  expect  that  some  flow  defect  will 
show  up  in  the  overall  calculations.  For  instance,  the  flow  area  may 
not  extend  radially  from  the  hub  completely  to  the  tip.  This  yields  a 
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flow  area  defect,  or  blockage,  which  is  theoretically  due  to  the  frictional 
effects  that  are  ignored  in  the  problem  but  may  also  include  measure¬ 
ment  errors.  Examples  of  the  frictional  effects  include  the  boundary 
layer  blockage  at  the  trailing  edge,  annulus  wall  friction,  and  the  fric¬ 
tion  at  the  planes  of  entrance  to  and  exit  from  the  control  volume.  The 
flow  defect  is  a  result  of  all  these  factors,  and  each  component  is  in¬ 
distinguishable  from  the  other.  Therefore,  the  blockage  factor  is  char¬ 
acteristic  of  the  friction  in  the  general  flow  process  and  not  necessarily 
representative  of  the  actual  boundary  layer  blockage  at  the  trailing  edge. 

From  consideration  of  Eqs.  (19  ),  (24'),  (27"),  (31),  and  (33  )  and 
Eqs.  (59),  (60),  and  (61)  describing  the  functional  form  of  the  variables, 
the  following  quantities  must  be  specified  by  solution  of  the  equations  in 
each  streamtube:  r^,  ^q,  bi,  b2,  gi,  g2»  dj.,  d2-  Since  there  are 
five  equations,  only  five  of  these  quantities  can  be  determined.  The 
following  section  outlines  the  selection  of  the  proper  five  constants  to  be 
determined  and  the  method  of  solution. 


2.7  SELECTION  OF  UNKNOWNS  AND  METHOD  OF  SOLUTION 

Substitution  of  Eq.  (61)  into  Eq.  (27  *)  leads  to  the  following  quad¬ 
ratic  formulation  of  Eq.  (27")  in  terms  of  d]_: 

d  1  V  t  -  2d  1  d2\'2  t-  ^2  \  3  -  mis  (27")  =0  (62) 


where 


R1IS  (27")  -  D'  + 


(r  -II  - 


r40) 


(r.tj 


'40^ 


4  “  *3 


k  —  — 

(v  277 


'30 


r31 


'30 

[>3  '  tli  dr  -  j  p3  th  dr 
'31 


r3  0 


f 


P3  '<•' 


'31 


,3  3  ■ 

~  I'3J  ~  r30 


2  ,  f}4  ,2 

- ('3|  -  r30) - ('31 '41 

7  4  ~  3  7  4  “  A  3 


2 

r30'40)  - 


-  ('3  ['41 


2 

'30'40 


k; 


)  - - - (r3]'4|  -  '30rt0^  - 

'  i  “  '-3  7i  ~  A3 


and  where 


'30  / 

Lvi,  V2.  V3l  J  1*3  oos-  y3  coS2  ^1  - 

'31 


n  lh  \ 
2rrr  / 


Ir,  r2,  r3l  dr 
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The  rearrangement  of  Eq.  (27")  into  Eq.  (62)  is  very  important  be¬ 
cause,  as  will  be  shown,  it  eliminates  the  necessity  of  iteration  on  one 
of  the  five  equations  since  dj  can  be  determined  exactly,  satisfying 
Eq,  (62),  if  the  other  variables  are  known  or  guessed. 

To  begin  the  solution,  the  particular  streamtube  is  selected  for 
which  the  overall- relative  total  pressure  loss  was  a  minimum.  The 
variables,  $3,  SU173,  M3',  are  assumed  constant  over  the  entire 

radial  height  of  this  initial  streamtube,  i.  e.  ,  b2  =  g2  =  ^2  =  0.  Thus, 
for  the  initial  streamtube,  d^  may  be  determined  from 

H]  =  .  IRHS  (27 ")/V \\A  (62') 

and  there  remains  exactly  five  undetermined  constants,  rgQ,  r3j_,  bj., 
gj,  di,  to  be  found  by  solution  of  the  five  governing  equations, 

Eqs.  (19  ),  (24  ),  (63'),  (31),  and  (33  ).  The  iteration  process  leading 
to  solution  is  illustrated  in  Fig.  7,  and  the  convergence  criterion  guar¬ 
antees  at  least  six  significant  digits  for  each  variable. 

The  solution  of  the  initial  streamtube  provides  boundaries  of  that 
streamtube  and  the  values  which  £3,  sin  73,  and  M3'  must  have  at  these 

boundaries  in  the  neighboring  streamtubes.  It  is  demanded  that  the 
solution  values  be  continuous  across  the  common  boundary  between  two 
streamtubes.  Then,  the  following  conditions  must  be  satisfied  for  every 
streamtube  except  the  initial  one: 


t>2  =  “  l>|Vrn 

(63) 

g2  -  tsin  >13  ~  g|)/7B 

(64) 

d2  =  (M'b  ~  d  1  )/r  13 

(65) 

where  rg,  j3g,  sin  7g,  and  Mg  are  determined  from  the  solution  in  a 
neighboring  streamtube.  Since  one  of  the  streamtube  bounding  radii  is 
determined  from  the  solution  of  a  neighboring  streamtube,  the  number 
of  unknown  quantities  has  been  reduced  from  eight  to  seven.  The  three 
additional  relations,  Eqs.  (63),  (64),  and  (65),  give  a  total  of  eight  equa¬ 
tions  so  that  the  problem  is  not  solvable  unless  one  of  the  five  governing 
equations  is  dropped.  The  axial  component  of  the  moment  of  momentum, 
Eq.  (33’),  was  eliminated  since  it  really  supplies  no  more  information 
to  the  system  of  equations  than  the  tangential  component  of  the  momentum 
equation,  Eq.  (24'). 
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Guess  r^Q 


<Do  current  guesses  of  all  variables  satisfy 
radial  component  of  the  momentum  equation? 

Convergence  Criteria:  1  x  10~5 
^  Yes 

<Do  current  guesses  of  all  variables  satisfy  the 
tangential  component  of  the  momentum  equation? 

Convergence  Criteria:  1  x  10” 7 
|  Yes 

<Do  current  guesses  of  all  variables  satisfy  the 
continuity  equation? 

Convergence  Criteria:  1  x  10_o 
|  Yes 

<Do  current  guesses  of  all  variables  satisfy  the 
moment  of  momentum  equation? 

Convergence  Criteria:  1  x  10-® 

Ives  . . . . 

I  Solution  obtained  in  the  initial  streamtube  I 


Fig.  7  Iteration  Scheme  for  the  Initial  Streamtube 
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If  Eq.  (65)  is  substituted  into  Eq.  (62),  the  following  relation  results 

A  dj  +  B  dj  +  C  =  0 


where 


.  0  V2  ,  Vs 
A  ■=  Vi  -  2  -  3  — r 

rD  rB 


B  =  2Mb 


,  r v 2  va 

? - 5“ 

L  rB  rB  j 


C  = 


Mb  2  — L  -  Ri  15(27  ") 


with  the  solution 


di 


-B  ±  vH*  -  4AC 
2A 


(62  ) 


Experience  has  shown  that,  in  general,  the  positive  sign  on  the  radical 
is  applicable  when  the  upper  boundary  of  the  streamtube  is  known  and 
the  negative  sign  applies  when  the  lower  boundary  is  known. 


The  iteration  cycle  yielding  solution  of  Eqs.  (19  ),  (24  ),  (62”), 
and  (31)  is  illustrated  in  Fig.  8.  The  convergence  factors  guarantee  at 
least  six  significant  figures  for  each  variable. 


It  remains  to  comment  here  that,  even  though  the  scheme  is  straight 
forward,  simple,  and  relatively  certain  of  finding  a  solution  if  it  exists, 
the  method  is  time  consuming.  One  complete  case  takes  approximately 
15  min  of  computation  time  for  the  IBM  360/50  computer  in  accomplish¬ 
ing  7000  or  more  iteration  steps  to  obtain  reasonable  accuracy.  There¬ 
fore,  a  much  faster  scheme  must  be  sought  if  the  problem  is  to  be 
economically  solved.  The  present  iteration  scheme  serves  only  to 
demonstrate  that  the  problem  can  be  solved  and  to  acquire  a  few 
example  results. 


2.8  DETERMINATION  OF  THE  WEIGHTS  FOR  EVALUATION  OF  THE  VOLUME  INTEGRALS 

Obtaining  and  correlating  values  for  the  weights  Wj(i  =  1,  2,  3,  4) 
in  Eqs.  (50),  (52),  (54),  and  (56)  presents  a  formidable  problem  in 
itself.  In  essence,  the  task  is  to  find  values  of  Wj  which,  when  substi¬ 
tuted  into  Eqs.  (50),  (52),  (54),  and  (56),  give  good  approximations  of 
the  values  of  R'i,  R'2,  ^  3>  and  R  4  obtained  by  more  accurate  evalua¬ 
tions  of  Eqs.  (22),  (23),  (26),  and  (35).  These  evaluations  may  involve 
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g2  =  (sin  Yb  - 


Solve  axial  component  of  the  momentum  equation  for 
J2  -  <*^-  d!>/rB 

/  Do  current  values  of  all  variables  satisfy  the  \ 
radial  component  of  the  momentum  equation? 

—  5 

\ _  Convergence  Criteria:  1  x  10 _  J 


Do  current  values  of  all  variables  satisfy  the 
tangential  component  of  the  momentum  equation? 

_  7 

Convergence  Criteria:  1  x  10 


Do  current  values  of  all  variables  satisfy  the 
continuity  equation? 

5 

Convergence  Criteria:  1  x  10 


Solution  obtained  in  any  streamtube  except  the  initial, 


Fig.  8  Iteration  Scheme  for  Solution  of  All  Streamtubes  except  the  Initial  Streamtube 
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measurements  at  other  planes  z,  zg  <  z  <  Z4.  However,  in  fluctuating 
flow  fields,  such  as  those  experienced  behind  a  rotor,  measurement 
interpretation  becomes  questionable.  Therefore,  even  the  more  accu¬ 
rate  evaluations  of  the  volume  integrals  probably  only  represent  low 
order  approximations. 


With  this  in  mind  consider  the  following  simple  approach.  For  the 
particular  case  under  consideration,  measurements  of  the  flow  proper¬ 
ties  are  available  2  in.  downstream  of  the  rotor  trailing  edge.  This  is 
interpreted  as  the  plane  z  =  Z4.  In  a  few  cases,  measurements  are 
simultaneously  taken  at  a  distance  of  1/2  in.  behind  the  rotor.  This 
plane  will  be  specified  z  =  ze.  Assuming  a  straight-line  variation  of 
the  integrals  between  the  z  stations  allows  the  approximation  of  the 
volume  integrals  in  the  following  manner  for  axisymmetric  flow: 


Hi  a  = 


S—  f '\rn\ 


V  K  R  I  \‘  T  3 


J  P3M3  cos  ya  sin /3a  ^1  +  — — -  rdr 

'31 


cos  ye 


5.- 

— -  Xg  <Z4  -  23)  H  -  zP)j| 


K>'0 


where 


rc0  /  .  \;> 

Xp  =  J  peMa  sin  fie  (l  ,  M;2j  rdr 

re  I 


and  where  X'  is  given  by  Eq.  (51) 


«. .  sl  i  «^.  r 


Rsa 


2kB 


1 


p .(i-^  *l) (-£)'* 


y;  (z3  -  <■  Y'f/e 


where 


(66) 


(67) 


(68) 


/  1  \ 

Ye  -  f  PeA  ■  ~ -  M'2Jr2dr 

r  is  I  ' 

and  where  Y*  is  given  by  Eq.  (53) 

|  ~r  J  P3M3  s'n  +  — - - M'a2^  -  - - ^  rdr 


(69) 


«3A  - 


VkR 


3  '31 


+ 


1  r .  a, 

— sin  ye  -  Ze 

VT'  A  e 

1  I-  c  geo 


(z  |  -  Z.\)  +  2'  UR 


(70) 
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where 


r,'°  /  \ 
-  /  ppm;  (l  f  m,,'2) 

re  I 


rdr 


(71) 


and  where  Z'  is  given  by  Eq.  (55) 

(  '«-*■  130 
r,!l 


n{,\ 


_ _  i 

VM»  / 


~ — -  J"  P3M3  sin  yj  +  — -  M 3 —  - — — ^  r 2d r 


1 


,r  A 


C  ««'■ 


Ze'U3 


z.^i  +  Z"  (zP 


(72) 


where 


J  Pc'Mr  ^  Me2^  '  r"dr 


(73) 


and  where  Z"  is  given  by  Eq.  (57).  In  each  case  T'  represents  an 
average  of  Tg  and  T4  which  should  partially  account  for  heat  transfer. 


A  certain  inaccuracy  is  admitted  in  these  equations  over  and  above 
the  measurement  uncertainty.  ‘  This  report  considers  tests  for  which 
the  pitch  angle  has  not  been  measured.  If  the  plane  z  ?  is  sufficiently 
far  downstream,  it  is  expected  that  the  direction  of  the  flow  would  be  in 
the  average  direction  of  the  bounding  walls.  For  the  case  of  cylindrical 
bounding  walls,  74  =  0,  however,  at  the  plane  z  =  ze,  such  an  approxi¬ 
mation  may  not  be  realistic.  Instead,  the  pitch  angle  (ye)  is  assumed 
constant,  hence  its  removal  from  the  integrals  of  Xg  in  Eq.  (66),  Zg  in 
Eq.  (70),  and  Z”  in  Eq.  (72),  and  it  is  to  be  related  to  the  average 
direction  of  the  streamtube  walls  at  the  plane  z  =  ze. 


By  using  Fig.  6,  the  following  relations  may  be  derived: 

«  tip  cos  yc  c  os  Pr 

M.'r 


tan  )> 


sin  ye  •=■ 


K. 

m;, 

"m; 


where  ye  is  the  angle  between  Mg  and  the  projection  of  Mg  into  a  tangen¬ 
tial  plane  and  where  7ez  is  the  angle  between  MgZ  and  the  projection  of 
Mg  into  a  radial  plane.  These  relations  lead  to 
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tan  ye  = 

tan  yez  cos 

tan  yC7.  cos  /3e 

V'l  +  tan2  yez  cos2  fie 

1 

v'l  +  tan2  yez  cos2  £e 

For  the  purposes  of  this  approximation,  j3e  in  these  relations  is  taken 
as  the  mass-averaged  relative  flow  angle  in  each  streamtube,  and  7ez 
is  given  by 

1  re 0  ~  r30  rel  “  r31  r40  “  re0 

tan  ye7.  =  — - +  -  +  -  + 

4  ze  —  zq  zc  -  z0 

Unfortunately,  the  addition  of  the  estimated  value  of  7e  results  in 
only  about  1 -percent  correction  of  the  continuity  balance  so  that  the  re¬ 
maining  error  may  be  due  to  inaccuracies  in  flow  property  measure¬ 
ments.  It  is  now  the  problem  to  force  the  continuity  equation  to  balance 
through  the  planes  ze  and  Z4,  but  it  is  not  feasible  to  apply  theoretical 
correction  factors  to  each  of  the  flow  property  measurements.  There¬ 
fore,  it  is  reasoned  that  the  unsteady  flow  at  the  plane  z  =  ze  may  be 
represented  as  the  periodic  passing  of  conditions  of  zero  mass  flow  and 
the  full  mass  flow  per  blade  passage.  The  sum  total  of  the  mass  flow 
per  passage  is  set  equal  to  the  total  mass  flow  measured  at  the  plane 
z  =  Z4.  This  procedure  will  then  yield  a  blockage  factor  (Ae/Aegeo)  to 
be  applied  to  the  integrals  evaluated  at  z  =  ze.  Essentially,  the  assump¬ 
tion  is  that  the  blocked  area  of  zero  mass  flow  is  so  small  in  terms  of 
time  and  the  time  constants  for  the  sensors  are  so  long  that  the  effects  of 
the  zero  mass  flow  sections  is  not  noted  in  the  measurements. 

Certain  errors  may  appear  with  the  use  of  this  correction  factor. 
These  are  related  to  the  fact  that  7  has  been  ignored  in  the  tests  to  date. 
Both  the  total  mass  flow  and  the  radial  locations  of  the  stream  surfaces 
at  z  =  ze  could  be  affected.  The  effect  should  be  small  since  the  mass 
flow  depends  on  the  cosine  of  this  angle. 

It  is  noted  that  the  correction  factor  directly  results  from  the  conti¬ 
nuity  equation  and  that  it  is  used  on  certain  terms  of  momentum  conser¬ 
vation  equations.  Although  an  analogous  correction  factor  could  be 
derived  from  the  one -dimensional  momentum  equation,  such  a  factor 
accounts  for  shear  forces  and  may  be  highly  sensitive  to  the  shearing 
effects  that  naturally  occur  in  this  flow  problem.  The  only  real  solution 
to  this  problem  is  the  correction  of  every  measurement  to  account  for 
the  very  turbulent,  unsteady  flows  that  must  occur.  Such  a  task  is 
beyond  the  scope  of  this  report. 


35 


AEDC-TR-71  -1 1  1 


Since  R'ia  has  been  obtained  by  Eq.  (66),  it  is  a  simple  matter  to 
set  this  value  equal  to  Eq.  (50)  and  then  to  solve  for  Wj_.  The  weights 
W2*  W3,  and  W4  may  be  found  in  a  similar  manner.  This  process  has 
guaranteed  that  the  volume  integrals  are  equal;  however,  it  has  not  en¬ 
sured  the  equality  of  the  solutions  by  both  methods.  This  fact  is  brought 
about  by  the  new  approximations  of  the  terms  B'  of  Eq.  (21)  and  E '  of 
Eq.  (29)  which  have  to  be  made  since  r3Q»  reQ,  an^  r40>  as  well  as  r3i, 
rej,  and  r4j,  may  not  be  coline ar  points.  However,  since  the  difference 
between  the  static  pressure  on  hub  and  tip  walls  is  not  very  large,  the 
differences  in  the  solutions  should  be  small.  At  any  rate,  it  is  time- 
wise  more  advantageous  to  directly  calculate  the  weights  by  equality  of 
volume  integrals  rather  than  through  an  iteration  process  for  the  weights 
guaranteeing  equality  of  solutions. 

In  general,  it  is  unlikely  that  the  representation  of  the  volume 
integrals  given  by  Eqs.  (66),  (68),  (70),  and  (72)  represent  more 
than  a  first-order  evaluation.  Improvement  of  this  scheme  can  come, 
first,  through  actual  measurement  of  the  pitch  angles,  and  second, 
through  providing  more  measurements  at  various  planes  between  Z3 
and  Z4.  Certainly,  more  elaborate  representations  of  the  volume  inte¬ 
grals  may  be  made  with  the  available  data.  For  instance,  a  quadradic 
approximation  of  the  integrands  may  be  hypothesized  rather  than  the 
linear  variation  assumed  here;  however,  the  basic  uncertainty  in  the 
interpretation  of  the  measured  data  lends  doubt  that  the  more  complex 
approximation  would  increase  accuracy. 


SECTION  III 

RESULTS  OF  CALCULATION 


The  following  paragraphs  present  the  calculation  results  for  rotor 
R1C2  to  exemplify  some  of  the  data  which  can  be  obtained  from  solution 
of  the  two-dimensional  sudden  area  expansion  flow  process.  The  rotor 
is  described  in  Bef.  2,  and  the  measured  data  are  summarized  in  that 
report.  Because  of  the  excessive  computation  time  required,  it  was 
necessary  to  limit  the  investigation  to  one  rotor  and,  in  general,  to  the 
cases  of  maximum  pressure  ratio  operating  conditions  for  each  speed 
line.  These  conditions  are  the  most  interesting  because  maximum 
efficiency  generally  occurs  at  these  conditions  for  the  rotors  tested  thus 
far.  One  speed  line  has  been  selected,  however,  for  which  calculations 
were  performed  for  all  operating  conditions  in  order  to  examine  the 
variation  of  sudden  area  expansion  loss  with  absolute  total  pressure 
ratio.  A  description  of  the  quasi-two-dimensional  solution  is  provided 
in  Appendix  I  for  comparison  purposes. 
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3.1  DETERMINATION  OF  THE  WEIGHT  NUMBERS 

The  weight  numbers  are  determined  such  that,  in  each  streamtube, 
the  volume  integrals  calculated  using  a  weighted  average  of  the  surface 
integrals  at  z  =  Z3  and  z  =  Z4  are  equivalent  to  the  volume  integrals 
calculated  by  assuming  straight-line  variation  of  the  surface  integrals 
with  axial  position  between  the  rotor  trailing-edge  plane  and  the  plane 
of  extra  measurements  and  between  the  plane  of  extra  measurements 
and  the  downstream  measuring  plane.  This  procedure  is  possible  only 
when  the  necessary  extra  flow  property  measurements  are  available; 
and,  at  present,  this  corresponds  to  the  maximum -pressure-ratio 
operating  conditions  only. 

Figure  9  is  a  plot  of  the  calculated  weight  numbers  in  function  of 
the  radial  position  of  the  streamtube  centers  at  the  downstream  meas¬ 
uring  plane.  Results  have  shown  that  W3  and  W4  are  very  nearly  equiv¬ 
alent.  This  could  have  been  expected  since  both  determine  essentially 
the  same  parameter;  i.  e. ,  the  tangential  component  of  the  volume  inte¬ 
gral  for  the  force  due  to  Coriolis  acceleration  (see  Eqs.  (54)  to  (57)). 

First-order  dependence  of  the  weight  numbers  on  the  proximity  of 
the  hub  and  tip  walls  is  evident.  This  probably  results  from  the  influ¬ 
ence  of  annulus  wall  friction  on  the  measured  data.  It  is  noted  that  W3 
and  W4  are  least  influenced  by  wall  proximity. 

The  weight  numbers,  especially  Wj_  and  W2,  are  noted  to  be  widely 
scattered.  They  exhibit  no  apparent  dependence  on  such  parameters  as 
wheel  speed,  Mach  number  at  the  downstream  measuring  station,  etc. 

It  was  originally  intended  to  correlate  these  weight  numbers  in  order 
to  attempt  calculations  at  conditions  for  which  the  extra  flow  property 
measurements  do  not  exist.  The  correlation  could  be  valuable  for 
calculations  at  maximum  pressure  ratio  operating  conditions  when  it  is 
desired  to  interpolate  between  constant  speed  characteristics,  but  it 
may  not  be  applicable  to  constant  speed  characteristics  when  the  abso¬ 
lute  total  pressure  ratio  is  not  maximum.  A  correlation  applicable  to 
the  latter  case  would  be  of  more  significance  for  this  purpose,  but  it 
is  not  possible  to  guarantee  the  reliability  of  such  a  correlation  using 
the  data  currently  available.  Therefore,  the  correlation  has  not  been 
attempted. 

It  has  been  pointed  out  that  there  will  .probably  exist  a  difference 
between  the  solutions  obtained  when  the  extra  plane  of  measurements 
are  used  in  calculating  the  volume  integrals  and  when  the  weight  num¬ 
bers  are  used.  These  differences  were  suggested  to  be  due  mainly  to 
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A  100%N//e  O  70%N/ZF 

□  90%N/ZiT  u  60%N/ZiT 

O  80%N/Z^  V  50%N/Ze 


Fig.  9  Variation  of  Weight  Numbers  as  a  Function  of  the  Radial  Position 
of  the  Streamtube  Centers  at  the  Downstream  Measuring  Plane 


differences  in  pressure- induced  forces  resulting  from  the  two  geometri¬ 
cal  approximations  of  the  streamtube  bounding  walls.  Figure  10  com¬ 
pares  the  results  of  solution  using  the  extra  plane  of  measurements  with 
the  results  using  the  calculated  weight  numbers  for  the  90-percent  N/\6 
data  point.  This  plot  is  typical  of  such  results,  and  it  is  noted  that 
major  differences  in  the  two  solutions  occur  in  the  streamtube  nearest 
the  tip  where  the  geometrical  differences  in  the  streamtubes  are  the 
greatest.  As  the  hub  wall  is  approached,  the  solutions  become  more 
nearly  the  same. 
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In  the  presentation  of  the  general  results  which  follows,  thie  results 
derived  from  the  solution  are  examined  using  the  weights  rather  than  the 
solution  determining  the  weights  because  only  this  alternative  is  avail¬ 
able  when  the  extra  plane  of  measurements  does  not  exist.  It  is  appar¬ 
ent  from  Fig.  10  that  this  solution  produces  results  that  are  less  erratic 
in  the  streamtube  nearest  the  tip. 


-  Solution  determining  weight  numbers 

-  Solution  using  weight  numbers 


20  30  40  50  60 


Fig.  10  Comporison  of  Solutions  Determining  Weight  Numbers  with  Solutions 
Using  Weight  Numbers  at  90-percent  N/\/0 
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3.2  GENERAL  RESULTS  OF  TH  E  TWO-DIMENSIONAL  SOLUTION 

Solution  of  the  two-dimensional  sudden  area  expansion  problem 
yields  the  radial  variation  of  the  relative  Mach  number  (Mg),  the  rela¬ 
tive  yaw  angle  (£3),  and  the  relative  pitch  angle  (73)  as  well  as  an  esti¬ 
mate  of  the  stream  surface  radial  locations  at  the  trailing-edge  plane. 
Figure  11  presents  the  radial  variation  of  M3  and  £3*  the  deviation 
angle  63  (i.  e.  ,  the  difference  between  £3  and  the  blade  angle  at  the 
trailing  edge),  and  73  for  various  rotational  speeds  of  rotor  R1C2  at 
maximum  absolute  total  pressure  ratio  operating  conditions.  These 
results  should  represent  circumferentially  averaged  values  at  each 
radial  station  in  the  trailing-edge  plane.  The  relative  Mach  number  is 
seen  to  decrease  with  decreasing  rotational  speed,  and  this  is  consistent 
with  the  decrease  in  the  relative  Mach  number  at  the  wheel  inlet. 

The  relative  yaw  angle  is  more  easily  understood  through  its  rela¬ 
tion  to  the  exit  blade  angle.  The  deviation  angle  ( 6 )  depicts  such  a 
representation,  and  it  should  be  positive  since  conditions  of  overturning 
by  blading  are  rarely  observed  in  practice.  Within  the  limitations  of 
the  problem  concept,  this  practical  condition  is  very  nearly  met  at  all 
rotational  speeds  except  100-percent  N /yfd.  Evidently,  the  flow  field 
at  design  speed  is  such  that  the  basic  assumptions  of  the  problem  are 
grossly  violated. 

The  results  for  the  relative  pitch  angle  indicate  that  this  angle  in¬ 
creases  with  decreasing  wheel  speed.  Such  a  result  is  extremely  sur¬ 
prising  since  the  radial  forces  in  the  wheel  are.  not  expected  to  be  suf¬ 
ficiently  high  to  justify  this  result.  It  will  be  shown  in  a  later  section 
that  this  discrepancy  can  be  directly  related  to  neglecting  frictional 
effects. 

The  results  for  the  radial  shifting  of  streamtubes  indicate,  in 
general,  that  all  streamtube  centers,  except  the  one  nearest  the  tip, 
are  located  between  their  known  positions  upstream  and  downstream  of 
the  rotor,  but  at  higher  rotational  speeds  {70-percent  N/\i0  and  above), 
most  of  the  shift  apparently  occurs  within  the  rotor.  Again  in  reference 
to  the  higher  speeds,  it  appears  that  the  streamtube  center  nearest  the 
rotor  tip  is  even  further  depressed  than  is  indicated  by  the  streamtube 
center  locations  at  the  measuring  plane  2.  0  in.  downstream  of  the  rotor. 
The  magnitude  of  the  depression  that  must  occur  within  the  rotor  de¬ 
creases  with  wheel  speed.  These  results  can  be  directly  related  to  the 
use  of  the  extra  plane  of  measurements  located  0.  5  in.  downstream  of 
the  rotor,  but  they  also  infer  that  the  difficulty  in  the  flow  field  indi¬ 
cated  by  the  wall  static  pressure  measurements  downstream  of  the  rotor 
(see  Refs.  2  or  9  or  Fig.  26  of  Ref.  8)  probably  results  from  separated 
flow  w-ithin  the  rotor. 
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Fig.  11  Solution  Results  lor  Rotor  R1C2  lor  Various  Rotor  Speeds 
at  Maximum  Absolute  Total  Pressure  Ratio  Operating 
Conditions 


The  main  objective  of  these  calculations  has  been  the  determination 
of  the  loss  in  relative  total  pressure  attributable  to  the  mixing  process 
which  occurs  when  the  flow  leaves  the  rotor  blade  passages.  The  calcu¬ 
lated  flow  properties  at  the  rotor  trailing  edge  allow  the  estimation  of 
an  average  relative  total  pressure  in  each  streamtube,  and  this  average 
value  is  used  to  calculate  the  loss  by  the  equation: 


The  radial  variation  of  sudden  area  expansion  loss  is  shown  in  Fig.  12. 
There  is  a  direct  correspondence  between  wgE  *n  this  figure  and  the 
solution  for  sin  73  shown  in  Fig.  11  demonstrating  that  Wgg  is  extremely 
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dependent  on  the  radial  component  of  the  flow  velocity.  However,  it  is 
very  doubtful  that  the  radial  component  behaves  as  the  present  solutions 
indicate.  This  point  will  be  examined  in  a  later  section. 

A  lOO%N//0  O  8O%N//0  U  6O%N//0 

□  90%N//6  O  70%N//6  V  50%N//6 


Fig.  12  Calculated  Radial  Variation  of  for  Rota  R1C2  for  Various 
Wheel  Speeds  at  Maximum  Absolute  Total  Pressure  Ratio 
Operating  Conditions 

/ 

Figure  13  shows  the  variation  of  w$e  as  a  one -dimensional  param¬ 
eter  for  conditions  of  maximum  absolute  total  pressure  ratio  in  com¬ 
parison  to  the  variation  of  the  overall  relative  total  pressure  loss. 

Three  points  are  shown  for  each  rotational  speed.  These  indicate  the 
repeatability  of  the  measurements  since  the  three  points  were  taken  at 
approximately  the  same  operating  conditions  and  emphasize  the  small 
effect  of  repeatability  on  the  solutions. 

The  results  shown  in  this  figure  indicate  that  wgjg  is  grossly  over¬ 
estimated  in  the  lower  speed  range  since  ugE  cannot  physically  be 
greater  than  u '.  This  fact  and  the  direct  relationship  between  ugE 
and  73  tend  to  indicate  that  the  solutions  for  73  are  physically  unreal¬ 
istic  in  the  lower  speed  ranges. 

It  is,  of  course,  un certain  that  any  of  the  two-dimensional  solutions 
presented  here  are  valid  approximations  to  the  real  flow  process.  Since 
the  overall  relative  total  pressure  loss  is  composed  of  sudden  area 
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u' 

Fig.  13  Comparison  of  the  Annulus  Average  Relative  Total  Pressure  Loss  Coefficient 
due  to  Sudden  Area  Expansion  with  the  Measured  Overall  Loss  at  Moximum 
Absolute  Total  Pressure  Ratio  Operating  Conditions 

expansion  loss,  passage  shock  loss,  and  profile  loss,  it  is  certain  that 
all  solutions  below  80-percent  N /  y[0  are  invalid  because  in  these  cases 
the  sum  of  the  calculated  sudden  area  expansion  loss  and  a  reasonable 
estimate  of  the  shock  losses  occuring  within  the  blade  passage  is  greater 
than  the  measured  overall  loss.  Even  the  80 -percent  N /\[6  case  is  very 
close  to  violating  this  condition.  However,  in  an  effort  to  determine  the 
influence  of  the  absolute  total  pressure  ratio  for  constant  wheel  speed 
operation,  the  solution  for  the  90 -percent  N fyf 9  speed  line  has  been  ob¬ 
tained.  To  perform  the  calculations,  the  weight  numbers,  determinable 
only  at  the  maximum  total  pressure  ratio,  were  assumed  constant  for 
the  entire  speed  line.  Example  solution  results  are  shown  in  Fig.  14, 
and  the  general  behavior  of  the  solutions  appears  quite  reasonable. 
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Figure  15a  shows  the  radial  variation  of  the  sudden  area  expansion, 
loss  for  the  same  cases  as  shown  in  Fig.  14,  and  in  Fig.  15b  the  one- 
dimensional  average  losses  are  shown  for  the  complete  speed  line  in 
comparison  with  the  overall  loss.  This  figure  shows  that  the  radial 
variation  of  use  changes  considerably  in  function  of  absolute  total  pres¬ 
sure  ratio,  whereas  the  one -dimensional  average  of  use  remains  about 
constant.  The  fact  that  the  average  is  about  constant  may  depend  on  the 
assumption  of  constant  weight  numbers.  Quasi-two-dimensional  solu¬ 
tions  for  this  problem  suggest  increasing  sudden  area  expansion  with 
decreasing  absolute  total  pressure  ratio  (see  Appendix  I). 

3.3  DISCREPANCIES  IN  THE  RESULTS  AND  THEIR  POSSIBLE  CAUSES 

Several  discrepancies  have  been  noted  to  result  from  solution  of 
the  two-dimensional  sudden  area  expansion  problem.  The  major  dis¬ 
crepancy  noted  in  the  present  solutions  is  that  sudden  area  expansion 
losses  become  extremely  large  at  the  lower  rotor  speeds.  This 
appears  to  be  directly  linked  to  the  rapid  increase  noted  in  the  pitch 
angle  73  as  rotor  speed  decreases.  It  was  further  noted  that  73  in¬ 
creases  in  the  vicinity  of  the  casing  walls  and  that  73  is  positive  while 
the  stream  surface  boundaries  generally  have  negative  slope  at  the 
trailing-edge  plane.  Such  effects  may  occur  at  discrete  points  in  the 
flow  field,  but  these  were  not  expected  to  be  representative  of  the  mean 
flow. 

The  result  of  increasing  73  with  decreasing  rotor  speed  may  be 
shown  to  result  from  the  particular  relation  between  the  radial  and 
axial  momentum  changes.  The  present  solution  is  a  result  of  the 
assumption  that  the  fluid  particles,  in  traversing  the  distance  between 
the  blade  trailing  edge  and  the  downstream  measuring  station,  are 
affected  only  by  the  radial  forces  due  to  pressure,  to  centripetal  accel¬ 
eration,  and  to  Coriolis  acceleration.  Figure  16a  shows  the  magnitude 
of  each- of  these  effects  for  the  entire  control  volume  at  each  rotational 
speed  as  well  as  the  net  radial  force  acting  on  the  control  volume.  This 
illustrates  that  in  the  present  problem  the  radially  inward  forces  due  to 
pressure  and  Coriolis  acceleration  are  larger  than  the  outward  force  due 
to  centripetal  acceleration.  Thus,  if  the  fluid  particles  are  to  leave  the 
control  volume  on  cylindrical  surfaces  (i.  e. ,  74  =  0),  then  the  radial 
momentum  at  the  trailing-edge  plane  must  be  directed  outward. 

The  axial  force  on  the  fluid  particles  in  the  control  volume  is  simply 
a  result  of  the  axial  static  pressure  differential.  The  ratio  of  the  resul¬ 
tant  radial  force  to  the  axial  force  is  shown  in  Fig.  16b.  The  ratio  of 
momentum  changes  are  equivalent.  Thus,  as  rotor  speed  decreases. 
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a.  Variation  of  the  Radial  Force  on  the  Control  Volume  as  a  Function 
of  Rotor  Speed 


b.  Relotion  of  the  Resultant  Radial  Force  to  the  Axial  Farce  on  the 
Control  Volume  as  a  Function  of  Rotor  Speed 

Fig.  16  Variotion  of  Control  Volume  Forces  as  a  Function  of  Rotor  Speed 
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the  radial  momentum  becomes  increasingly  larger  than  the  value  of  the 
axial  momentum.  Therefore,  the  average  radial  Mach  number  in  the 
control  volume  is  predicted  larger  in  value  than  the  average  axial  Mach 
number,  and  73  will  increase  correspondingly. 

At  this  point,  it  is  reasonable  to  question  the  applicability  of  the 
equations  in  the  absolute  coordinate  system  to  this  problem.  The  change 
in  the  absolute  radial  momentum  depends  only  on  the  force  due  to  the 
annulus  pressure  differential  under  the  assumptions  of  steady,  friction¬ 
less  flow  of  a  weightless  gas.  Figure  16a  shows  that  the  force  due  to 
pressure  is  greater  than  the  resultant  force  in  the  relative  coordinate 
system.  Thus,  the  radial  momentum  change  must  be  even  greater  in 
an  absolute  coordinate  system  than  the  present  solution  indicated  for  the 
relative  coordinate  system.  Therefore,  the  discrepancies  would  be  even 
greater. 

The  present  solution  neglects  the  forces  due  to  friction,  unsteady 
flow,  and  the  weight  of  the  gas.  The  detailed  consideration  of  friction 
and  unsteady  flows  is  beyond  the  scope  of  this  report,  whereas  the  weight 
of  the  gas  may  still  be  considered  negligible.  However,  it  must  be  ques¬ 
tioned  whether  any  of  these  could  produce  the  necessary  correcting 
forces. 

Of  the  three  neglected  forces,  the  force  due  to  friction  is  the  only 
one  considered  here  since  it  is  probably  the  largest.  Although  solutions 
have  been  obtained  under  the  assumption  of  frictionless  flow,  the  meas¬ 
urements  used  to  obtain  the  solution  result  from  a  real,  viscous  process. 
Thus,  the  present  solutions  suggest  no  more  than  the  distribution  of  the 
flow  properties  at  the  trailing  edge  if  the  measured  flow  properties  re¬ 
sulted  from  an  inviscid  process. 

Friction  in  the  relative  coordinate  system  will  cause  fluid  particles 
near  the  casings  to  tend  to  adhere  to  the  casing  walls  which  are  moving 
in  this  system.  This  causes  the  relative  velocity  to  approach  the  velocity 
of  the  walls;  i.  e. ,  a  velocity  directed  opposite  to  the  wheel  rotation  with 
the  magnitude  of  the  wheel  speed.  This  result  has  two  related  conse¬ 
quences  on  the  present  solution,  indicating  that  friction  in  the  measure¬ 
ments  is  retained  by  use  of  the  inviscid  equations.  First,  since  the 
radial  component  of  the  force  due  to  Coriolis  acceleration  is  a  function 
of  the  tangential  component  of  the  relative  velocity,  this  radial  force, 
directed  inward,  increases  in  magnitude  very  rapidly  as  the  wall  is 
approached.  The  solutions  for  73  reflect  this  phenomenon  since  73 
tends  to  increase  as  either  casing  wall  is  approached. 
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Second,  the  relative  Mach  number  near  the  tip  casing  wall  should 
increase  in  value,  ultimately  becoming  equal  to  the  Mach  number  of  the 
wheel.  This  condition  seems  to  be  occurring  in  the  solutions  for  the 
trailing-edge  plane. 

In  the  vicinity  of  the  tip  casing  wall,  the  addition  of  the  shear  stress 
term  to  the  radial  component  of  the  momentum  equation  will  tend  to 
counteract  the  large  increase  in  the  force  due  to  Coriolis  acceleration 
thereby  reducing  the  net  resultant  force  inward.  In  the  vicinity  of  the 
hub  wall,  the  local  velocity  gradients  with  respect  to  r  are  opposite  to 
those  at  the  tip  and,  thus,  tend  to  increase  the  net  inward  force.  How¬ 
ever,  if  real  flows  are  to  be  considered  at  the  trailing -edge  plane,  it  is 
necessary  to  include  the  boundary  condition  that  the  relative  velocity 
tends  to  zero  in  the  blade  passage  as  the  hub  is  approached.  This  sig¬ 
nificantly  alters  the  velocity  gradient  at  the  trailing-edge  plane  and  re¬ 
quires  that  the  radial  component  of  the  force  due  to  Coriolis  accelera¬ 
tion  approach  zero.  Thus,  the  net  resultant  force  inward  should  tend  to 
decrease  at  the  hub. 

It  has  been  qualitatively  shown  that  the  effect  of  friction  will  tend  to 
reduce  the  discrepancies  noted  in  the  solution  for  the  radial  component  of 
the  relative  Mach  number  in  the  vicinity  of  the  casing  walls.  This  im¬ 
plies  73  may  be  expected  to  behave  more  rationally  if  the  flow  of  a  real 
fluid  is  considered.  However,  it  is  obvious  in  Fig.  11  that  the  discrep¬ 
ancies  are  continued  into  the  main-stream  flow. 

In  Appendix  II,  an  attempt  to  estimate  the  effects  of  free  turbulent 
flows  has  been  made  to  show  that  free  turbulent  shear  cannot  be  ne¬ 
glected  on  a  stream  surface.  The  analysis  shows  that  the  effect  of  fric¬ 
tion  on  the  radial  component  of  the  momentum  equation  is  not  sufficient 
to  reduce  the  discrepancies  in  the  center  portion  of  the  annulus.  It 
demonstrates  that  the  tangential  component  of  the  momentum  equation  is 
primarily  affected  by  free  turbulence.  The  force  due  to  Coriolis  accel¬ 
eration  then  produces  the  coupling  necessary  to  actually  correct  the 
present  solutions  for  the  radial  component  of  velocity. 

The  analysis  in  Appendix  II  demonstrates  that  it  is  necessary  to  con¬ 
sider  the  fully  three-dimensional  nature  of  the  problem.  The  tangential 
variation  of  velocity  in  the  blade  passage  initialize  irregular  velocity 
variations  in  the  blade  wakes  which  result  in  unbalanced  shear  forces  in 
the  wakes.  Radial  secondary  flows  in  these  wakes  are  shown  to  be  ex-  ■ 
tremely  important.  The  analysis  is  only  an  approximation,  however, 
because  of  the  severe  limit  in  knowledge  about  such  flow  phenomenon. 
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Neglecting  friction  in  compressible  flows  is  commensurate  with 
neglecting  heat  transfer  in  the  energy  equation.  However,  it  must  be 
pointed  out  that  heat  transfer  apparently  does  occur  from  the  fluid  to 
the  tip  casing.  This  statement  is  based  on  the  fact  that  measurements 
at  the  extra  measuring. plane  0.  5  in.  behind  the  rotor  and  the  measur¬ 
ing  plane  2.  0  in.  behind  the  rotor  indicate  that  the  entropy  of  the  fluid 
decreases  between  these  planes  in  the  streamtube  nearest  the  tip  for 
all  rotational  speeds.  Sufficient  temperature  measurements  allowing 
calculation  of  the  heat  transfer  rate  have  not  been  made  in  tests  to  date. 

Because  of  the  great  number  of  assumptions  incorporated  into  this 
analysis,  it  is  impossible  to  point  out  any  other  single  feature,  besides 
friction,  which  has  had  a  large  influence  on  the  present  results.  The 
author  has  investigated  the  consequences  of  errors  in  the  weight  num¬ 
bers  and  in  the  specification  of  the  static  pressure  at  the  hub  in  the 
trailing-edge  plane.  Such  errors  are  relatively  minor  in  comparison 
to  the  gross  error  introduced  by  neglecting  friction.  Certain  other 
assumptions,  such  as  straight-line  static  pressure  variation  radially 
across  the  annulus,  the  conical  form  of  the  stream  surfaces,  and  others 
which  purport  to  estimate  unmeasured  flow  properties,  may  be  deemed 
critical  when  it  is  possible  to  include  friction  in  the  present  analysis. 

The  interpretation  of  the  existing  measurements  may  ultimately  be  ques¬ 
tioned.  It  should  be  recalled,  however,  that  frictional  effects  were  de¬ 
picted  as  large  in  Appendix  II  under  the  basic  assumption  that  radial 
velocity  components  were  relatively  large  in  the  blade  wakes;  i.  e. ,  con¬ 
sideration  of  secondary  flows  becomes  an  essential  part  of  the  solution. 
Furthermore,  it  is  the  variation  of  flow  properties  from  blade  to  blade 
that  initiates  the  velocity  gradients  of  sufficient  magnitude  to  produce  the 
large  frictional  effects.  Thus,  it  can  be  concluded  that  the  sudden  area 
expansion  problem  for  the  blunt  trailing-edge  rotor  is  a  fully  three- 
dimensional  problem  which  will  not  be  completely  understood  without 
a  general  and  complete  three-dimensional  analysis. 


SECTION  IV 

CONCLUDING  REMARKS 


The  sudden  area  expansion  process  is  fully  three-dimensional  in 
rotor  flows.  Tangential  blade -to -blade  flow  property  variations  and 
radial  velocities  in  the  blade  boundary  layers  and  wakes,  indicative  of 
secondary  flows,  play  an  extremely  important  part  in  determining  the 
flow  development.  At  present,  these  are  not  well  enough  understood  to 
allow  any  more  than  a  qualitative  estimate  of  their  effects. 
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This  report  has  investigated  the  sudden  area  expansion  process  as 
a  quasi -three -dimensional,  frictionless,  adiabatic  flow  process  by  the 
use  of  the  streamtube  approximation,  and  it  has  been  demonstrated  that 
the  real  flow  process  cannot  be  approximated  in  such  a  manner.  The 
general  behavior  of  the  solutions  for  the  relative  Mach  number  and  the 
relative  yaw  angle  at  the  trailing  edge  may  be  considered  quite  good; 
however,  the  behavior  of  the  solutions  for  the  relative  pitch  angle  must 
be  considered  very  poor  as  a  description  of  the  average  flow  process. 

The  reasons  for  the  difficulties  are  relatable  to  neglecting  friction. 
Friction  at  solid  boundaries  is  important,  but  it  is  shown  that  the  pri¬ 
mary  effect  of  friction  on  a  problem  approached  by  the  streamtube 
approximation  is  relatable  to  free  turbulent  shear.  Velocity  profiles 
are  highly  nonaxisymmetric  at  the  trailing -edge  plane  because  of  the 
blade  element  curvature.  This  results  in  unbalanced  shear  forces  in 
the  blade  wakes.  The  solutions  of  the  problem  considered  in  this  report 
indicate  that  these  effects  must  be  accounted  for.  In  addition,  sizable 
radial  velocities  in  the  wakes  are  also  required  to  account  for  the  dis¬ 
crepancies  noted  in  the  present  solutions. 

It  is  believed  that  the  discrepancies  occurred  because  real-flow 
measurements  were  used  with  the  ideal-flow  equations.  Measurements 
were  available  to  describe  the  flow  properties  downstream  of  the  rotor, 
and  the  flow  conditions  at  the  rotor  trailing  edge  were  demanded  from 
the  inviscid  flow  equations.  The  resulting  solutions  can  only  be  inter¬ 
preted  as  the  conditions  at  the  trailing  edge  if  the  process  leading  to  the 
measured  data  were  inviscid. 

The  quasi-two-dimensional  equations  lead  to  more  reasonable  re¬ 
sults  for  sudden  area  expansion  flow  process.  This  result  can  only  be 
related  to  the  real  flow  process  if  it  is  hypothesized  that  neglecting 
friction  and  radial  secondary  flows  are  at  least  partially  compensating 
assumptions.  It  should  be  emphasized  that  the  solutions  of  the  quasi- 
two-dimensional  equations  will  lead  to  incorrect  results  if  the  stream- 
tube  assumption  is  used. 

The  development  and  form  of  the  quasi-three -dimensional  equations 
show  that  it  is  not  currently  possible  to  solve  the  problem  for  the  flow 
conditions  downstream  of  the  rotor  if  only  the  conditions  at  the  rotor 
trailing  edge  are  supplied.  The  solution  requires  some  knowledge  of  the 
flow  process.  In  particular,  it  requires  the  static  pressure  develop¬ 
ment  along  the  annulus  walls  as  well  as  some  information  about  the  vol¬ 
ume  integrals  of  the  forces  due  to  Coriolis  and  centripetal  acceleration. 
Only  experience  with  the  solution  of  the  equations  in  a  form  similar  to 
the  present  analysis  will  lead  to  knowledge  necessary  to  solve  the  in¬ 
verted  problem. 
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An  experimental  and  theoretical  program  is  proposed  to  continue 
study  of  this  problem.  Some  suggested  experiments  are: 

1.  Sudden  area  expansion  in  pipe  flows,  similar  to  those  of 
Ref.  12,  to  include  a  sufficient  variety  of  area  ratios 
and  a  further  investigation  of  the  effect  of  geometry 
over  the  complete  range  of  initial  Mach  number 

(0  to  1.  0)  and  possibly  varying  back  pressure  when 
the  inlet  flow  is  choked  and  the  downstream  flow  be¬ 
comes  supersonic.  Such  an  investigation  is  aimed  at 
determining  the  static  pressure  on  the  trailing  sur¬ 
face.  The  theory  for  this  problem  is  simple  for  fric¬ 
tionless  flow  and  should  agree  very  well  with  experi¬ 
mental  results. 

2.  Flows  in  cascades  designed  to  produce  symmetrical 
outlet  flows  from  each  passage.  The  theory  is  still 
one -dimensional  and  the  effects  of  shear  stress  are 
probably  still  negligible,  but  the  objective  is  to  deter¬ 
mine  the  effects  of  wake  flows  on  the  trailing -surface 
static  pressure.  Studies  of  free  turbulence  may  be  of 
interest  to  determine  wake  development. 

3.  Flows  in  cascades  of  curved  surface  blading  produc¬ 
ing  tangential  gradients  in  the  initial  flow  before  the 
sudden  area  expansion.  The  theory  becomes  two- 
dimensional,  and  the  objective  is  to  determine  the 
applicability  of  basically  one -dimensional  theories 
such  as  Ref.  7.  Discrepancies  due  to  shear  stress 
are  expected  to  become  increasingly  important.  Studies 
of  free  turbulence  should  definitely  be  attempted. 

4.  Flows  through  curved  blades  in  annular  cascade  to 
determine  the  effect  of  radial  pressure  gradients. 

The  theory  must  become  completely  three-dimensional 
now,  and  very  significant  discrepancies  should  occur 
if  the  effect  of  shear  stress  is  neglected  in  the  analysis. 
Wake  studies  become  extremely  important  here  since 
this  is  the  simplest  case  through  which  experiments 
can  demonstrate  the  effects  of  radial  forces  on  the 
development  of  the  wakes. 

5.  Flows  in  rotors.  At  present,  because  of  the  lack  of 
techniques  for  measuring  flow  properties  in  the  rotating 
system,  this  study  must  be  completely  theoretical.  All 
the  information  gathered  by  the  previous  studies  may 
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still  be  insufficient  except  at  very  low  rotor  speeds 
where  the  forces  due  to  Coriolis  and  centripetal 
acceleration  may  be  considered  negligible.  This 
study  is  intimately  connected  with  a  study  of  the 
flow  development  through  the  blading  since  the  ad¬ 
vantage  of  knowing  the  initial  conditions  before  the 
sudden  area  expansion  process  does  not  exist  in 
the  case  of  rotor  flows. 
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APPENDIX  I 

THE  QUASI-TWO-DIMENSIONAL  APPROACH 


FORMULATION  OF  THE  PROBLEM  AND  REVIEW  OF  A  RELATED  EXPERIMENT 

In  an  effor  to  establish  the  approximate  order  of  magnitude  of  loss 
due  to  the  sudden  expansion  process,  the  analysis  in  Ref.  7  used  the 
momentum  theorem  with  certain  simplifying  assumptions.  The  main 
assumptions  were: 

1.  The  flow  is  steady,  adiabatic,  and  axisymmetric  in  the 
absolute  frame  of  reference  with  concentric  cylindrical 
stream  surfaces. 

2.  Friction  is  neglected. 

3.  Radial  secondary  flows  are  neglected. 

The  loss  in  total  pressure  then  becomes  a  function  only  of  the  area 
ratio  and  the  component  of  the  Mach  number  in  the  direction  of  the  area 
increase.  Thus  a  fully  three-dimensional  problem  is  reduced  to  solu¬ 
tion  in  one  dimension.  The  axial  component  of  the  momentum  equation, 
neglecting  the  shear  stress  term  and  the  weight  of  the  gas,  is  used  in 
the  form 


P4^z4A4  -  P3W23A3  =  p3A3  -  p4A4  +  p3'(A4-A3)  (I_l) 

where  pg'  represents  the  pressure  on  the  surface  of  the  blade  trailing 
edges.  In  Ref.  7,  P3'  is  assumed  equal  to  P3  which  leads  to 

p3(A4  -  A3)  +  p3A3  4-  P3W£3A3  -  P4A4  -  P4W24A4  =  0 


By  using  the  one -dimensional  continuity  equation,  the  energy  equation, 
and  the  equation  of  state  to  eliminate  P4/P3,  the  equation  above  becomes 


A4 

—  +.  K  M3  2  cos2  /3  3 
a3 _ 

1 3  cos  $3  ^  +  — 2 - 


I  -  K  M4 2  cos2 


(1-2) 


If  station  3  is  assumed  to  be  located  at  a  plane  such  that  the  flow  has 
not  yet  experienced  the  sudden  area  expansion  and  such  that  the  rotor 
may  not  further  impart  momentum  to  the  flow,  then  substitution  of  the 
two-dimensional  flow  triangle  relations  into  the  tangential  component 
of  the  momentum  equation  leads  to 
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M3  sin  /S3  -M4  sin  /34 


(1-3) 


since  the  stream  surfaces  are  assumed  cylindrical.  Thus  Eqs.  (1-2) 
and  (1-3)  relate  the  variables  A4/A3,  Mg,  ,83,  M'4,  and  34.  Any  three 
of  these  may  be  chosen  as  knowns  and  the  other  two  may  then  be  deter¬ 
mined  uniquely.  The  static  pressure  ratio  and  the  total  pressure  re¬ 
covery  are  then  determined  from  the  continuity  equation  in  the  form: 
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(1-5) 


The  one -dimensional  approximation  is  valid  as  long  as  no  cross - 
flow  components  exist;  however,  the  approximation  may  be  applied  to 
real  flows  if  the  velocity  profiles  are  fully  developed.  This  means  that 
integration  of  shear  stress  over  a  control  surface  perpendicular  to  the 
mean  flow  direction  must  be  zero;  then  neglecting  friction  merely  re¬ 
quires  that  friction  along  solid  boundaries  be  negligible. 

These  considerations  imply  that  the  quasi-two-dimensional  equations 
are  best  applied  to  the  entire  annulus  and  not  to  individual  streamtubes. 
The  results  for  an  extension  of  the  analysis  presented  in  the  main  body 
of  this  report  affirm  this  conclusion  and  point  out  that  tangential  blade - 
to-blade  velocity  variations  produce  large,  shear  forces  in  the  blade 
wakes  that  must  be  accounted  for  in  the  streamtube  approximation. 

Under  these  circumstances,  the  best  approximation  seems  to  be 
calculation  of  the  sudden  area  expansion  loss  based  on  average  flow 
properties  for  the  entire  annulus  and,  assuming  the  same  resulting  loss, 
occurs  in  each  streamtube. 

In  Ref.  12,  an  attempt  is  made  to  correlate  empirically  the  value  of 
P3'  in  Eq.  (I- 1)  with  the  actual  measured  losses  in  one -dimensional  pipe 
flow  (i.  e. ,  j3j  =  0,  j  =  3,  4)  for  both  sudden  area  expansion  and  contrac¬ 
tion.  The  assumptions  of  this  reference  and  those  of  Ref.  7  are  other¬ 
wise  the  same.  The  average  base  pressure  for  sudden  area  expansion 
and  contraction  is  measured  experimentally,  and  the  actual  losses  are 
compared  with  those  determined  analytically  using  the  average  measured 
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base  pressure.  A  good  correlation  results  indicating  that  use  of 
Eq.  (1-2)  yields  somewhat  lower  losses  than  those  which  occur  in  experi¬ 
ment.  If  p3'  is  left  in  Eq.  (1-1),  an  equation  similar  to  Eq.  (1-2)  results. 

~~  -  l)+  1  +kU32  c°s2^3  1  -  KM*'2  cos2/34 

M3  cos  /S3  (l  +  M32)  M4  cos  /34  (l  +  M  4  2) 

Since  Eq.  (1-3)  remains  unchanged,  the  Eqs.  (1-2')  and  (1-3)  relate  the 
variables  A4/A3,  M3,  jS3,  M4,  j34,  and  p3' /p3.  Thus,  it  is  necessary 
to  know  four  of  these  quantities  before  the  other  two  may  be  determined 
uniquely  However,  -if  sufficient  data  are  available  to  establish  an  accu¬ 
rate  correlation  of  P3'/P3  in  function  of,  say,  A4/A3  and  M3  cos  £3.  then 
again  selection  of  any  three  of  the  variables  as  known  quantities  allows 
unique  solution  for  the  remaining  two  variables. 

In  Ref.  12,  two  area  ratios  (outlet  area/ inlet  area)  were  tested, 

2.  25  and  4.  00,  using  both  circular  and  square  ducting.  Definite  depend¬ 
ence  of  base  pressure  on  geometry  was  noted.  The  results  for  expansion 
indicate  that,  as  the  Mach  number  before  expansion  approaches  sonic, 
as  much  as  10 -percent  error  may  be  made  in  the  calculation  of  the  total 
pressure  recovery  factor  through  the  use  of  an  equation  like  Eq.  (1-2) 
rather  than  Eq.  (1-2')  at  an  area  ratio  of  2.  25  and  about  15-percent 
error  occurs  in  the  case  of  an  area  ratio  of  4.  00.  The  difference  in  the 
results  of  Eqs.  (1-2)  and  (1-2')  decreases  as  Mach  number  decreases. 

These  experiments  indicate  that  some  error  is  made  in  the  use  of 
Eq.  (1-2);  i.  e. ,  by  the  assumption  that  P3'  =  P3  and  that  the  error  de¬ 
creases  as  both  the  Mach  number  before  expansion  and  the  area  ratio 
decrease.  Insufficient  data  exist  at  the  present  to  allow  correlation  of 
p3'  for  area  ratios  useful  in  this  investigation  (1.  37  to  1.  83). 

It  might  be  useful  to  make  additional  tests  extending  the  data  of 
Ref.  12  to  smaller  area  ratios  so  that  Eq.  (1-2')  may  be  used  in  the  one- 
dimensional  analysis  of  the  blunt  trailing-edge  blading.  The  reasons 
against  such  a  program  include  the  fact  that  the  flow  in  an  apparatus 
like  that  of  Ref.  12  may  be  quite  different  from  that  experienced  by  a 
blunt  trailing-edge  stator  .and  certainly  must  be  different  from  the  flow 
experienced  by  a  rotor  due  to  the  peculiarities  of  the  rotating  flow  field. 
However,  until  flow  measurement  in  a  rotating  flow  field  becomes  prac¬ 
tical,  researchers  must  be  satisfied  with  experiments  performed  in  a 
nonrotating  field  conscientiously  applied  as  guides  for  describing  the 
flow  in  a  rotating  coordinate  system.  Certainly  refined  experiments 
similar  to  that  in  Ref.  12  should  be  attempted  for  flow  in  annular 
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cascades  when  the  problem  is  critically  important  to  the  analysis,  but 
for  the  moment,  the  simplest  experiment  should  be  attempted  to  deter¬ 
mine  the  order  of  magnitude  of  the  error  made  in  using  Eq.  (1-2)  rather 
thanEq.  (1-2'). 

METHODS  OF  SOLUTION  AND  RESULTS  OF  CALCULATION 

The  particular  use  of  EqS.  (1-2)  to  (1-5)  for  the  calculations  of  this 
report  requires  some  further  rearrangement  of  those  equations.  Meas¬ 
urements  at  the  downstream  measuring  plane  may  be  used  to  yield  aver¬ 
age  values  of  M4  and  184  for  the  annulus.  Thus,  one  of  the  remaining 
three  variables,  A4/A2,  M3,  or  £3,  must  be  specified  in  order  to  deter¬ 
mine  a  solution  of  Eqs.  (1-2)  and  (1-3).  It  is  immediately  obvious  that 
M3  cannot  be  specified  and  that  A4/A3  should  not  be  specified  since  one 
of  the  objectives  of  the  calculation  was  to  determine  the  effective  flow 
area  at  the  trailing  edge.  Under  the  basic  .assumption  of  the  applicability 
of  these  equations,  it  is  necessary  that  the  flow  should  not  be  separated 
at  either  the  tr ailing-edge  plane  or  the  downstream  measuring  plane. 
Therefore,  the  flow  angle  (£3)  should  be  reasonably  close  to  the  exit 
blade  angle  (83).  Their  equivalence  is  assumed  for  the  calculations  for 
Method  One  below. 

A  second  possible  method  of  calculation  requires  somewhat  more 
leniency  in  both  data  interpretation  and  the  assumptions.  This  method 
involves  the  specification  of  the  static  pressure  ratio  accomplished  by 
the  process.  Measurements  of  static  pressure  are  made  along  the  hub 
and  tip  casing  walls  between  the  rotor  trailing- edge  plane  and  down¬ 
stream  measuring  stations  except  at  the  hub  wall  in  the  trailing-edge 
plane  where  measurement  is  not  physically  possible.  If  this  static  pres¬ 
sure  were  available  together  with  the  static  pressure  measurement  at 
the  rotor  tip,  then  a  linear  variation  could  be  used  to  approximate  the 
static  pressure  at  any  radius.  Use  is  made  of  the  Legrange  interpolation 
formula  and  all  static  pressures  measured  on  the  hub  wall  to  extrapolate 
for  the  hub  wall  static  pressure  at  the  rotor  trailing-edge  plane.  The 
closest  measured  pressure  is  0.  25  in.  downstream  of  the  rotor  trailing 
edge,  and  this  is  one  of  five  used  in  the  extrapolation,  the  most  distant 
being  at  the  plane  of  the  downstream  flow  measurements.  This  is 
Method  Two  described  below. 


Description  and  Results  of  Method  One 


The  critical  Mach  number,  given  by 

M"2  - 


!1Z_LM2 


K  -  1 


M2  +  I 
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may  be  used  to  put  Eq.  (1-3)  into  the  form 

M3*  sin  £3  -  M4  *  sin  /3 4 

which,  for  given  j3g,  jS^.,  and  M4,  may  be  solved  for  M3*  leading  to 
M3.  Equation  (1-2)  may  then  be  solved  for  A4/A3.  These  computations 
yield  sufficient  information  to  determine  the  static  and  total  pressure 
ratios  from  Eqs.  (1-4)  and  (1-5).  The  relative  total  pressure  loss  attrib¬ 
utable  to  the  sudden  area  expansion  process  is  then  given  by 


The  average  one-dimensional  losses  are  plotted  versus  absolute  total 
pressure  ratio  in  Fig.  1-1  for  rotor  R1C2. 

Careful  examination  of  Fig.  1-1  shows  that  sudden  expansion  loss 
generally  increases  as  the  absolute  total  pressure  ratio  decreases  for 
constant  speed  operation  except  in  the  case  of  design  speed  operation. 
Since  maximum  pressure  ratio  implies  maximum  diffusion  to  a  minimum 
passage  exit  Mach  number,  most  of  the  results  show  reasonable  varia¬ 
tion  for  the  loss  due  to  the  sudden  area  expansion  flow  process.  As  the 
pressure  ratio  is  decreased  at  constant  wheel  speed,  the  trailing-edge 
Mach  number  increases,  implying  that  the  level  of  sudden  expansion 
loss  should  increase. 

However,  each  of  the  curves  for  design  speed  operation  exhibit  an 
abrupt  decrease  in  loss  at  minimum  total  pressure  ratio.  This  result 
of  the  calculations  can  be  directly  linked  to  the  occurrence  of  a  sharp  in¬ 
crease  in  the  axial  Mach  number  downstream  of  the  rotor  when  back 
pressure  is  reduced  to  a  minimum.  The  calculated  trailing-edge  rela¬ 
tive  Mach  numbers  tend  to  show  that  the  flow  through  the  rotor  has 
actually  been  accelerated  for  these  minimum-pressure-ratio  operating 
conditions.  Static  pressure  measurements  at  the  rotor  tip  also  tend  to 
confirm  the  acceleration  at  minimum  pressure  ratio  (see  Refs.  2  and  3 
and  Fig.  26  of  Ref.  8). 


Description  and  Results  of  Method  Two 


Equation  (1-4)  may  be  solved  for  A4/A3  and  used  to  eliminate  the 
area  ratio  from  Eq.  (1-2).  Then  substitution  of  Eq.  (1-3)  leads  to 

kM  '2  sin  04  cos  £4 
tanp3  =  - - - 

1  —  kM  2  cos  2  /3 4 

p  4 
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from  which  £3  may  be  calculated  when  P3/P4,  M4,  and  £4  are  given. 

The  calculation  procedure  is  then  similar  to  Method  One  leading  to  the 
calculation  of  the  relative  total  pressure  ratio  by  Eq.  (1-5).  The  aver¬ 
age  one -dimensional  sudden  expansion  losses  are  plotted  versus  abso¬ 
lute  total  pressure  ratio  in  Fig.  1-2  for  rotor  R1C2. 

The  general  trends  of  variation  of  sudden  expansion  loss  for  each 
speed  line  in  Fig.  1-2  is  very  similar  to  that  observed  in  Fig.  1-1;  i.  e.  , 
increasing  loss  with  decreasing  pressure  ratio.  The  calculation  results 
must  then  be  considered  to  yield  plausible  results.  However,  it  is  noted 
that  the  magnitude  of  sudden  expansion  loss  is  considerably  reduced  with 
the  use  of  Method  Two. 

Figure  1-3  shows  the  variation  of  the  deviation  angle  (63)  at  the 
trailing-edge  plane  in  function  of  absolute  total  pressure  ratio  for 
selected  wheel  speeds.  The  deviation  angle  is  defined  as  the  difference 
between  the  flow  direction  and  the  blade  angle,  both  measured  in  a  plane 
perpendicular  to  the  radial  direction.  The  deviation  angle  is  a  measure 
of  the  flow  guidance  accomplished  by  the  blading.  In  Method  One,  flow 
guidance  was  assumed  perfect;  i.  e.  ,  63  was  identically  zero. 


COMPARISON  OF  METHODS 

In  this  section,  it  is  intended  to  compare  and  contrast  the  results  of 
the  three  methods  (two  quasi-two-dimensional  methods  and  the  quasi- 
three -dimensional  method)  to  determine  which,  if  any,  can  be  used  to 
approximate  the  losses  due  to  the  sudden  area  expansion  process  that 
occurs  as  the  flow  leaves  the  rotor  blade  passages. 

From  the  discussion  of  Fig.  13,  it  must  be  concluded  that  sudden 
area  expansion  loss  is  grossly  over-estimated  by  the  quasi-three- 
dimensional  calculations  at  low  rotor  speed;  however,  the  possibility  of 
reasonable  predictions  at  high  speeds  must  be  considered.  The  behavior 
of  the  solutions  at  90-percent  in  function  of  absolute  total  pressure 

ratio  shown  in  Figs.  14  and  15  may  be  reasonable  even  though  the  trend 
of  loss  variation  is  not  substantiated  by  the  quasi-two-dimensional 
results.  The  single  feature  which  must  lend  doubt  to  the  accuracy  of  the 
solutions  is  the  large  values  of  73  and  the  fact  that,  in  most  cases,  the 
flow  at  the  blade  trailing  edge  is  not  directed  par  ailed  to  the  bounding 
stream  surfaces.  This  occurrence  has  been  directly  linked  to  neglect¬ 
ing  viscous  shear.  (It  might  be  commented  here  that  an  attempt  has  been 
made  to  force  solution  of  the  equations  to  negative  values  of  73.  Such 
solutions  were  not  possible,  indicating  that  the  inviscid  equations  could 
not  be  solved  if  73  were  approximated  by  the  local  slope  of  the  stream 
surfaces  at  the  trailing-edge  plane. ) 
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Open  symbols  -  traverse  measurements 
Closed  symbols  -  rake  measurements 

A  100%N//6  O  80%N//9  P  60%N//6 

D  90%N//6  O  70%N//F  V  50%N//9 


Fig.  1—2  Sudden  Area  Expansion  Loss  Calculated  by  Method  Two  versus 
Absolute  Total  Pressure  Ratio 
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Open  symbols  -  traverse  measurements 
Closed  symbols  -  rake  measurements 


A  lOO%N//0 
□  9O%N//0 


O  8O%N//0 
O  7O%N//0 


V  60%N//e 

V  50%N//e 


1-3  Variations  of  the  Deviation  Angle  Calculated  by  Method  Two 
as  a  Function  of  Absolute  Total  Pressure  Ratio 
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Since  much  of  the  effort  to  correct  the  quasi-three -dimensional 
solutions  in  the  previous  section  was  directed  toward  reduction  of  the 
radial  velocity,  it  appears  reasonable  that  the  quasi-two-dimensional 
solutions  which  ignore  the  existence  of  a  radial  component  in  velocity 
should  give  more  reasonable  results.  However,  neglecting  shear  stress 
would  be  just  as  critical  in  the  case  of  the  quasi-two-dimensional  solu¬ 
tions  as  it  is  in  the  case  of  the  quasi-three-dimensional  solutions  if  the 
streamtube  approximation  were  used. 

An  approximation  similar  to  the  quasi-two-dimensional  method  may 
be  attempted  using  the  quasi-three-dimensional  equations.  This  is 
accomplished  by  solving  the  entire  annulus  as  if  it  were  the  initial  stream- 
tube  (see  Section  2.  7).  Essentially,  this  means  all  solution  variables 
are  assumed  constant  at  the  trailing-edge  plane,  but  that  both  radial 
and  tangential  components  of  velocity  are  allowed  to  occur.  The  results 
for  sudden  area  expansion  loss  are  shown  in  comparison  to  the  com¬ 
pletely  quasi-three -dimensional  results,  the  quasi-two-dimensional  re¬ 
sults  by  both  methods,  and  the  measured  overall  loss  in  Fig.  1-4.  This 
method  produces  a  significant  reduction  in  the  estimated  loss  but  not 
nearly  sufficient  to  consider  it  a  correct  approximation  to  the  actual  loss 
at  low  rotor  speeds. 

Comparison  of  the  two  quasi-two -dimensional  methods  reveals  that 
Method  One  generally  yields  higher  relative  total  pressure  losses  for  the 
process  and  smaller  values  for  the  flow  angle,  static  pressure,  and  flow 
area  at  the  trailing-edge  plane.  The  tendancy  of  increasing  sudden  ex¬ 
pansion  loss  with  decreasing  absolute  total  pressure  ratio  for  constant 
speed  operation  is  apparent  in  both  quasi-two-dimensional  methods. 

This  is  more  consistent  with  ideas  of  how  flow  in  the  passage  behaves 
than  the  results  of  the  quasi-three -dimensional  analysis  which  suggests 
that  the  loss  is  constant. 

Considering  the  relation  between  the  two  quasi-two-dimensional 
methods  it  may  be  concluded  that  the  results  would  be  the  same  by  each 
method  if: 

1.  The  flow  were  assumed  to  have  a  small,  generally  posi¬ 
tive  deviation  from  the  blade  angle  at  the  trailing  edge 
for  use  in  Method  One,  or 

2.  The  extrapolation  procedures  necessary  for  Method  Two 
were  to  have  produced  a  hub  wall  static  pressure  at  the 
trailing  edge  that  was  too  large,  or  both  of  these. 
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Fig.  1—4  Comparison  of  the  Measured  Overall  Relative  Total  Pressure  Loss  vvith 
the  Results  of  vario.us  Methods  Predicting  Sudden  Area  Expansion  Loss 
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Either  of  these  alternatives  could  have  occurred  for  any  given  set 
of  measured  conditions  downstream  of  the  rotor,  thereby  making  the 
choice  between  them  quite  arbitrary.  Ultimately,  the  decision  as  to 
which  method  is  more  realistic  for  this  problem  can  be  made  only  in 
connection  with  the  flow  model  as  a  whole.  This  report  does  not  deal 
with  the  consequences  of  selecting  the  calculation  procedure;  however, 
some  comments  on  these  two  alternatives  may  be  made  here. 

First,  some  deviation  from  perfect  flow  guidance  certainly  must 
be  expected  at  the  trailing  edge.  Its  magnitude  could  be  calculated 
from  Method  Two  if  the  trailing-edge  static  pressure  were  known  with 
some  certainty  and,  of  course,  if  the  assumptions  of  the  analysis  were 
obeyed  in  reality.  Second,  it  is  felt  that  the  extrapolated  static  pres¬ 
sure  on  the  hub  wall  is  not  predicted  too  large.  As  a  matter  of  fact, 
in  a  small  number  of  cases  (actually  amounting  to  only  1  percent  of  the 
total  number),  the  extrapolated  pressure  was  too  small  to  allow  solu¬ 
tion  of  the  equations. 

This  qualitative  analysis  suggests  that,  if  the  assumptions  of  the 
theory  are  obeyed.  Method  Two  is  probably  superior  to  Method  One.  If 
the  assumptions  are  not  obeyed,  in  particular,  if  the  flow  contains  radial 
flow  components,  then  the  results  of  the  quasi-three -dimensional  analy¬ 
sis  suggests  that  Method  One  is  more  reasonable  since  it  predicts 
larger  losses  in  relative  total  pressure. 
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APPENDIX  II 

AN  ESTIMATE  OF  THE  EFFECT  OF  FREE  TURBULENT 
FLOWS  ON  THE  STREAMTUBE  APPROXIMATION 


According  to  Ref.  13  and  others,  turbulent  boundary  layer  theory 
assumes  that  the  shear  stress  may  be  given  by 

/  \  <*T 

T  =  pW  +  ()  — 

oy 

where  v  is  the  laminar  kinematic  viscosity,  e  is  the  virtual  kinematic 
viscosity  of  turbulent  flows,  u  is  the  local  time-averaged  velocity  of  the 
mean  motion,  and  y  denotes  the  direction  transverse  to  the  mean  motion. 
For  free  turbulent  flows,  such  as  jets  and  wakes,  v  may  be  neglected 
in  comparison  to  e.  From  page  481  of  Ref.  13, 

.  .  .  L.  Prandtl  assumed  that  the  dimensions  of  the  lumps  of 
fluid  which  move  in  a  transverse  direction  during  turbulent 
mixing  are  of  the  same  order  of  magnitude  as  the  width  of 
the  mixing  zone.  .  .  .  The  virtual  kinematic  viscosity,  e,  is 
now  formed  by  multiplying  the  maximum  difference  in  the 
time -mean  flow  velocity  with  a  length  which  is  proportioned 
to  the  width,  b,  of  the  mixing  zone.  Thus 

f  =  M  b  (II- 1) 

Here  denotes  a  dimensionless  number  to  be  determined 
experimentally. 

The  value  of  e  is  often  assumed  to  remain  constant  over  the  whole 
width  of  the  wake  at  every  cross  section;  and,  in  the  case  of  plane  wakes. 
Ref.  14  demonstrates  that  e  is  constant,  independent  of  the  distance 
from  the  body,  as  long  as  velocity  profile  similarity  exists.  The  latter 
result  is  then  valid  only  at  large  distances  from  the  wake -producing  body, 
whereas  the  former  assumption  may  be  corrected  by  introducting  an 
intermittency  factor  accounting  for  the  periods  of  time  in  which  turbu¬ 
lence  occurs  at  a  particular  point  in  the  flow  field. 


Most  theoretical  wake  flow  studies  have  been  applied  to  incom¬ 
pressible  flows  with  negligible  pressure  gradients,  and  these  have  been 
successful  only  in  application  to  large  distances  from  the  body.  In 
application  to  the  present  problem,  Eq.  (II- 1)  for  e  and  the  resulting 
equation 

T  ~  P  *1  b  (umax  “  umin)  2) 
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for  shear  stress  will  be  used  to  arrive  at  estimates  for  these  quantities 
in  the  vicinity  of  the  rotor  trailing  edges.  The  quantity  (K^)  is  assumed 
a  constant,  but  in  reality  it  may  depend  on  the  type  of  problem  to  which 
the  approximation  for  e  is  being  applied  as  well  as  on  the  distance  from 
the  body  at  which  r  is  being  calculated. 

Certain  further  assumptions  must  be  made  for  the  present  applica¬ 
tion.  The  solutions  of  the  quasi-three -dimensional  sudden  area  expan¬ 
sion  problem  make  it  possible  to  estimate  the  average  relative  velocity 
and  density  at  the  plane  of  the  rotor  trailing  edge.  The  variation  of 
these  quantities  with  percent  of  design  rotor  speed  is  shown  in  Fig.  H-l 
(top).  By  assuming  that  the  wake  width  (b)  may  be  approximated  by  the 
blade  trailing-edge  thickness  (th)  and  that  the  velocity  difference  is  given 
by 

umax  —  umin  =  ^3  —  ® 


then 


(II-3) 


where  rref  is  the  radius  of  the  blading  reference  profile.  The  results 
of  this  calculation  are  shown  in  Fig.  H-l  (middle). 


Several  examples  using  the  approximation  for  e  expressed  by 
Eq.  (II- 1)  are  given  in  Ref.  13.  The  value  of  the  constant  apparently 
depends  on  the  problem.  The  following  results  have  been  noted: 

(0.014  frec-jet  boundary 
0.047  quadi-lhree-dimensional  wake 
-  ^0.0333  wake  behind  row  of  bars 

J0.0185  two-dimensional  jet 
f 0.0128  circular  jet 

According  to  Ref.  15,  may  be  on  the  order  of  0.  1  for  the  case  of  the 
mixing  of  parallel  streams  when  initial  velocity  gradients  on  the  body 
are  taken  into  account. 


For  the  purposes  of  this  investigation,  a  value  of  0.  0333  is  chosen 
for  because  the  wake  behind  a  row  of  bars  best  represents  the  prob¬ 
lem  and  because  it  is  convenient  to  state  that  the  final  results  may  be 
too  small  or  too  large  by  a  factor  of  three-.  Then,  e  is  given  by  the 
right-hand  scale  of  Fig.  H-l  (middle). 
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Two  surfaces  over  which  shear  stress  must  be  integrated  are  of 
particular  importance  to  the  present  problem.  The  first  consists  of  the 
planes  perpendicular  to  the  axial  direction.  For  such  a  plane  the  inte¬ 
grand  of  the  shear  stress  integral  in  the  momentum  equation  is  (-tz)rz 
which  may  be  resolved  into  radial  and  tangential  components: 

(-lz)  rz  =  (-ar)  rzr  +  (-afl)  rz0  (II -4) 


The  second  type  of  surface  is  the  conical  stream  surface  which  bounds 
each  streamtube.  By  resolving  the  force  which  produces  the  shear 
stress  on  this  type  of  surface  into  its  components,  it  is  found  that  (-t)T 
may  be  represented  by 


(-t)r  =  (-ar)  t./t  sin  £  +  (-a0)  (rr0  cos  £  +  tzq  sin  £)  h  (-^)  rrz  cos  £ 


(n-5) 


where  ?  is  the  inclination  of  the  stream  surface  with  respect  to  axial 
direction  measured  in  a  plane  containing  the  axis  and  where  the  normal 
stresses  rrr  and  rzz  have  been  neglected  in  the  transformation. 

By  using  Eqs.  (II-4)  and  (II— 5),  the  shear  stress  integral  of  the 
momentum  equation  may  be  written 

J  (-1)  i’dA  -  (~ar)  J  "zr  dA3  -  J*  ^zi-dA*  +  J"  ^zr  sin  £  dAff 

(A)  *-U3;  (A  4)  vAw> 


rz  0  <IA4 


'  ,/  TyQ  ^3  1  J  T 

MA3)  (A4) 

-  J  (rr0  cos  f  +  rz0sin£)dAw 

l  A  it  )  -J 


H  (-k)  J  rrzcosfdAw 


(II  -6> 


where  A3  is  the  surface  area  at  the  blade  trailing- edge  plane,  A4  is  the 
surface  area  at  the  downstream  measuring  plane,  and  A\y  represents 
the  area  of  the  axisymmetric  surfaces  bounding  the  streamtubes  between 
A3  and  A4,  and  where  from  page  54  of  Ref.  13, 


*7 


(II- 7a) 
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(II- 7b) 


(II- 7c) 


The  surface  A3  has  been  assumed  to  be  located  infinitesimally  far 
upstream  of  the  blade  trailing  edge.  Therefore,  the  integral  of  shear 
stress  over  this  surface  is  associated  with  the  boundary  layer  on  the 
blade  profiles.  Since  turbulence  is  inhibited  by  the  presence  of  solid 
boundaries,  shear  stress  is  expected  to  be  less  at  this  plane  than  that 
set  up  by  the  free  turbulence  occurring  downstream  of  the  rotor.  For 
this  reason,  the  integrals  over  A3  will  be  ignored  in  the  analysis  which 
follows. 


Consider  first  the  radial  component  of  Eq.  (II- 6).  The  surface 
Ayy  is 

r 4  tZi 

Aff  =  2^  J  rQ(z)  dz  +  2ir  J  rj(z)dz  =  2nA/.  (rQ  -  rj) 

*3  *3 

if  average  values  of  ro  and  ri  are  selected.  The  largest  possible  value 
of  sin  £  is  given  by 

sin  |f|  <  tan  'f|  =  - — 

az 


Then, 


A4  =  !7  (r$  -  if)  =  -j-  sin  |f|  Aw 


Thus,  if  an  average  value  of  rzr  is  selected  for  the  entire  control  volume, 

ar  ■  J  (-f)  rdA  =  ffr  =  -  rzr  A4  (1  ±  2) 

(A) 

where  the  sign  of  the  shear  stress  integral  over  Aw  depends  on  the 
slope  of  the  surfaces.  The  maxim  tom  value  of  this  integral  is  then 
-3rzrA4.  If  the  force  due  to  shear  stress  is  to  cancel  the  resultant 
radial  force  shown  in  Fig.  16a,  then  using  the  values  of  p  and  e 
(neglecting  v)  shown,  respectively,  in  Fig.  II- 1  top  and  middle  and 
Eq.  (II-7a)  for  Tzr,  it  is  possible  to  estimate  the  magnitude  of  the  term 
(9Wr/9z  +  9Wz/9r).  The  result,  as  shown  in  Fig.  n-1  (bottom),  indi¬ 
cates  that  this  term  should  be  approximately  2  x  10^  sec”*.  If  such 
gradients  occur,  then  Wr  will  approach  zero  so  that,  essentially,  9WZ/  9r 
alone  must  be  on  this  order  of  magnitude.  At  design  speed,  B  Wz/9r  is 
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on  the  order  of  104  in  the  center  of  the  annulus  at  the  downstream  meas¬ 
uring  plane,  and  it  becomes  much  less  than  this  as  rotor  speed  decreases. 
Thus,  gradients  on  this  order  (and  higher,  of  course^  are  observed  only 
in  the  proximity  of  the  casing  walls  so  some  other  influence  must  be 
critical. 

Consider  now  the  tangential  component  of  Eq.  (41)  involving  r rQ 
and  tzq  as  given  by  Eqs.  (II-7b)  and  (II-7c),  respectively.  In  the  center 
portion  of  the  annulus ,  9Wg/8r  may  be  ignored.  As  a  first  approxima¬ 
tion,  9W@/8z  may  be  ignored.  It  must  be  noted  that,  even  if  the  analysis 
demonstrates  that  Wr  becomes  small  in  the  free  stream,  neglecting  the 
change  of  Wr  in  the  tangential  direction  would  deny  the  existence  of 
secondary  radial  flows  which  may  occur  in  the  wakes.  Therefore,  this 
term  must-be  included  in  the  present  analysis,  and  the  tangential  com¬ 
ponent  of  Eq.  (II-6)  becomes,  approximately 

a  0  ’  J  (-t)rdA 

(A) 


=  Ue  »  -  f  Tzdd&4  -J  ( TjQ  cos  £  T  rz0  sin  dA  V; 
(A4)  (A*) 


•  >  1 

+  —  -  cos  c  +  —  -  sin  cl  dA\v 

r  dd  r  dd  / 


where  v  has  been  ignored  in  comparison  to  e.  The  term  Wg/r  is  ne¬ 
glected  in  comparison  to  the  velocity  change  in  the  tangential  direction. 
Since  e  is  of  large  magnitude  in  the  wakes  and  negligible  in  the  free 
stream,  this  equation  reduces  to 


p4e  (9w74  b 

fffl  = - - - — —  rref  —  (rQ  -  r|)  -  4rrrref 


*4  / 

J  H 

z  3 


pi  <9W, 

-  cos  c  +  - sin 

dd  dO 


Since 


then 


w2  =  w2  +  w l  W'2 


<?W  Wr  <9wr  ViQ  dViQ  W  z  dW, 


30  W  dd  W  dd  W  06 
or,  solving  for  8Wr/  80 

<9ffr  w  w,  <5wz 

ffr  dd  Wr 


•Of 


dz 


dd 


dd 


(II-8) 
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where  it  has  been  assumed  that  the  chance  of  W q  in  the  tangential  direc¬ 
tion  is  negligible  in  the  wake.  The  result  of  the  approximation  for 
3Wr/96  depends  on  the  local  value  of  Wr.  It  will  be  assumed  that,  in 
the  wake,  all  velocity  components  are  of  comparable  orders  of  magni¬ 
tude;  i.  e. , 

W2  =  3VC?  «  3\Vr2  •=  3W^ 

Then, 

to,  _  2  to, 

00  ~  00 

For  plane  wakes  at  large  distances  from  the  wake  producing  body. 
Ref.  14  demonstrates  that  b  is  proportioned  to  z  2  and  that  the  velocity 
defect  varies  as  z-'2,  so  that 

i  0\\  „ 

—  -  b  =  constant 

r  06 

If  b  may  be  represented  by  th,  if 

j_  __  ^3  ~  0 _ 2  w3 

r  do  J_  /th  \  ,  /  ih  \ 

2  rrcf  \  r~ )  ^'(t) 

and  if  Eq.  (II-3)  is  used  to  eliminate  e,  then  for  constant  density 

ff 0  =  -  477  P3K]  R!3rref  jh'0  -  T|)  +  2  (2  cos  £  +  sin  $  (z.t  -  z3)j 

-4!7p3K[W3rref  +  t  (2  cos  £  +  sin  £)1  (1  in.)  (II-9) 


since 

ro  -  rT  ~  1  in- 

z 4  -  7.3  =  2  in. 

Equation  (II-9)  shows  that  integration  of  the  shear  stress  over  the 
stream  surfaces  is  quite  important.  By  assuming  that  ?  is  small, 
commensurate  with  small  values  of  Wr  in  the  free  stream,  then  ff# 
becomes 

UO  =  -36  W3 rref  U  in-) 
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The  results  of  the  calculation  of  the  absolute  magnitude  of  the  force  due 
to  friction  in  the  tangential  direction  (ffg)  are  shown  in  Fig.  H-2  where 
the  computations  have  used  p 3  and  W3  given  in  Fig.  II- la, 

k  j  =  0.0333 
rref  =  20.450  in. 


%N//S‘ 

Fig.  11-2  Comparison  of  Forces  Affecting  the  Tongential  Component  of 
the  Momentum  Equation 

Consider  the  following  simplified  form  of  the  tangential  component 
of  the  momentum  equation  neglecting  unsteady  effects  and  the  weight  of 
the  gas: 

mW40  -  mW30  =  ff0  -  iConolis# 

or,  solving  for  W39  while  noting  W40  and  ffg  are  negative  in  the  present 
case, 

*30  =  ^fConolis^  -  »«|W4(9|  I-  |ff0|J  (II-  1 0) 

The  tangential  components  of  the  force  due  to  Coriolis  acceleration  and 
of  the  quasi-two-dimensional  average  momentum  passing  through  A4 
are  also  shown  in  Fig.  II-2  in  relation  to  this  equation.  It  is  seen  from 
Eq.  (II- 10)  that  the  addition  of  the  tangential  frictional  forces  will  cause 
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Wgg  to  be  less  negative.  This  is  indicated  in  Fig.  II- 2  by  the  differences 
between  the  tangential  momentum  at  A3  calculated  by  neglecting  friction 
and  by  the  inclusion  of  the  frictional  force.  Since  the  radial  component 
of  Coriolis  acceleration  depends  on  W q,  the  neglection  of  friction  implies 
an  overestimation  of  the  radial  force  due  to  Coriolis  acceleration. 

It  is  noted  that,  if  Wr  is  decreased  by  the  addition  of  friction,  then 
the  tangential  force  due  to  Coriolis  acceleration  is  also  reduced  and  that, 
if  W30  is  reduced,  it  is  possible  that  the  good  agreement  between  the 
present  solution  for  $3  and  the  blade  angle  may  be  destroyed.  Both  of 
these  factors  are  fundamentally  dependent  on  the  coupling  of  the  compon¬ 
ents  of  the  momentum  equation  through,  not  only  the  Coriolis  accelera¬ 
tion  term  which  connects  the  radial  and  tangential  components,  but  also 
the  introduction  of  friction  which  connects  all  three  components.  The 
effects  of  coupling  make  the  final  result  extremely  difficult  to  predict, 
but  it  is  felt  that  the  tendency  toward  correction  of  the  present  solutions 
has  been  satisfactorily  demonstrated. 
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